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Abstract 

ZFITTER is a Fortran program for the calculation of fermion pair production and 
radiative corrections at high energy e~^e~ colliders; it is also suitable for other 
applications where electroweak radiative corrections appear. ZFITTER is based on a 
semi-analytical approach to the calculation of radiative corrections in the Standard 
Model. We present a summary of new features of the ZFITTER program version 6.42 
compared to version 6.21. The most important additions are: (i) some higher-order 
QED corrections to fermion pair production, (ii) electroweak one-loop corrections 
to atomic parity violation, (iii) electroweak one-loop corrections to Pgt'e production, 
(iv) electroweak two-loop corrections to the W boson mass and the effective weak 
mixing angle. 
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New version summary 



Title of program: ZFITTER version 6.42 (18 May 2005) 

Autliors of original program: D. Bardin, P. Christova, M. Jack, L. Kali- 
novskaya, A. Olchevski, S. Riemann, T. Riemann 

Program obtainable from: 

littp://www-zeutlien.desy.de/tlieory/researcli/zfitter/ (main web site), 
/ af s / cern. ch / user /b /bardindy / public / ZF6_42 

Reference for ZFITTER version 6.21 D. Bardin et al., Comput. Phys. Com- 
mun. 133 (2001) 229-395 

Operating system: UNIX/LINUX, program tested under, e.g., HP-UX and 
PC/Linux 

Programming language used: FORTRAN 77 
High speed storage required: < 2 MB 

No. of cards in combined program and test deck: about 26,200 

Keywords: Quantum electrodynamics (QED), Standard Model, electroweak 
interactions, heavy boson Z, e"'"e~-annihilation, fermion pair production, ra- 
diative corrections, initial state radiation (ISR), final state radiation (FSR), 
QED interference, SLD, LEP, ILC. 

Nature of the physical problem: Fermion pair production is an important 
reaction for precision tests of the Standard Model, at LEP/SLC and future 
linear colliders at higher energies. For this purpose, QED, electroweak and 
QCD radiative corrections have to be calculated with high precision, includ- 
ing higher order effects. Multi parameter fits used to extract model parameters 
from experimental measurements require a program of sufficient flexibility and 
high calculational speed. ZFITTER combines these two aspects by employing 
analytical integrations of matrix elements and at most one one dimensional 
numerical integration, as well as a variety of flags deflning the physics content 
used. The calculated predictions are typically at the per mille precision level, 
sometimes better. 

Method of solution: Numerical integration of analytical formulae. 

Restrictions on the complexity of the problem: Fermion pair production is 
described below the top quark pair production threshold. Photonic corrections 
are taken into account with simple cuts on photon energy, or the energies and 
acoUinearity of the two fermions, and one fermion production angle. The treat- 
ment of Bhabha scattering is less advanced. 

Typical running time: On a Pentium IV PC installation (2.8 GHz) using 
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g77 under Linux 2.4.21, approximately 23 sec are needed to run the standard 
test of subroutine ZFTEST. This result is for a default/recommended setting 
of the input parameters, with all corrections in the Standard Model switched 
on. ZFTEST computes 12 cross-sections and cross-section asymmetries for 8 
energies with 5 interfaces, i.e. about 360 cross-sections in 23 seconds. 
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1 Introduction 

The Fortran program ZFITTER is based on a semi-analytical approach to the 
calculation of fermion pair production in e~^e~ annihilation at a wide range of 
centre-of-mass energies, including SLC/LEPl, LEP2, and ILC energies below 
the tt-threshold. ZFITTER allows the calculation of several quantities needed 
for precision studies of the Standard Model: 

• Mw - the W boson mass; 

• Tz, Tyy, - total (and also partial) Z and W boson decay widths; 

• da /d cos d - differential cross-sections; 

• (Jt - total cross- sect ions; 

• ApB ^ forward-backward asymmetries; 

• - left-right asymmetries; 

• Apou ^ final state polarisation effects for r leptons. 

All observables are calculated including radiative corrections using the (run- 
ning) fine structure constant a, the muon decay constant G^, the Z mass Mz 
as well as the fermion masses and the Higgs mass Mh as input. 

ZFITTER version 6.21 and higher was mainly intended for the use at LEPl/SLC 
energies, and versions since 6.30 are also adapted to the LEP2 kinematics. 
Various interfacing subroutines (short: interfaces) allow the user to calculate 
observables, e.g., for fits to the experimental data with different sets of free 
parameters. 

The Fortran package DIZET, a library for the calculation of electroweak radia- 
tive corrections, is part of the ZFITTER distribution. It can also be used in a 
stand-alone mode. On default, DIZET performs the following calculations: 

• by call of subroutine DIZET: W mass, Z and W partial and total decay 
widths; 

• by call of subroutine ROKANC: four weak neutral-current (NC) form factors, 
running electromagnetic and strong couplings needed for the calculation of 
effective NC Born cross sections for the production of massless fermions 
(however, the mass of the top quark appearing in the virtual state of loop 
diagrams for the process e^e~ ff is not ignored); 

• by call of subroutine RHOCC: the corresponding form factors and running 
strong coupling for the calculation of effective CC Born cross sections; 

• by call of subroutine ZU_APV: Q]y{Z,A) - the weak charge used for the 
description of parity violation in heavy atoms. 

If needed, the form factors of cross sections may be made to contain the contri- 
butions from WW and ZZ box diagrams thus ensuring the correct kinematic 
behaviour over a larger energy range compared to the Z pole. 
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ZFITTER version 6.21 [1] was released in July 1999 and was described in [2]. 
Since then, there were several important developments of the program. The 
current release is that of ZFITTER version 6.42, dated 18 May 2005. This article 
describes the changes and additions in ZFITTER from version 6.21 to version 
6.42. Other sources of information on ZFITTER and its use are given by the 
previous program description and the references therein, the ZFITTER webpage 
[3], the studies [4,5], the 1999/2000 CERN LEP2 workshop proceedings [6], 
and the studies of the LEPEWGG (LEP electroweak working group) [7]. 

The essential changes from ZFITTER v. 6.21 to v. 6.42 in terms of physics 
topics are: 

• Higher order QED corrections to fermion-pair production, of importance at 
energies off the Z boson peak; 

• Electroweak corrections to the weak charge Qw^ describing the parity vio- 
lation effects in atoms, of importance for so-called global Standard Model 
fits; 

• Electroweak corrections to Vf,Vf. production, of importance for a precise de- 
scription of vv^ production; 

• Electroweak two-loop corrections to M^y and the effective weak mixing angle 
gj^2,e//^^^ of importance for global Standard Model fits and for precise 
predictions of the Higgs mass Mn- 

Further, an option to change the strength of the Wtb vertex, \yth\i was imple- 
mented in DIZET since version 6.30: \yth\ is now one of the parameters of sub- 
routine DIZET. The default is the Standard Model calculation with \yth\ = 1- 
Several new interface routines deal with this case, see Appendix C for the 
details. 

This program description update is organized as follows. In Sections 2 to 5, we 
describe the improved or new physics issues. Section 6 contains a description 
of the input parameters and pseudo observables, and Section 7 reproduces 
reference outputs from a running of the sample program package. A Summary 
closes this update note. In Appendices A to C we collect, for the convenience 
of the user, updated ZFITTER and DIZET user guides and an overview of the 
user interfaces. The presentation assumes, of course, a familiarity of the reader 
with [2]. 

With the present update of ZFITTER and its description, the maintenance of 
the program has been migrated to a group of volunteers, the ZFITTER support 
group. The main webpage has been migrated to: 
http : // www-zeuthen. desy. de / theory / research / zfitter / 

The decision to do so was taken together with the original authors of the 
program. Since we did not want to change the list of authors in view of the 
long-standing history of the program, the idea to create a ZFITTER support 
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group was considered to be an appropriate way to handle the maintenance 
situation requiring permanent care about the program. The ZFITTER support 
group is composed of authors of ZFITTER version 6.21, long-term users of the 
program, and colleagues who contributed substantially to its present state. We 
hope that a kernel of us will stay with ZFITTER as long as the code is needed 
by the community. 

A correct citation of the ZFITTER package will include reference [2], together 
with the present update. 



2 Higher order QED corrections 



Measurements at LEPl/SLC were performed in the vicinity of the Z resonance 
so that many of the photonic corrections were suppressed. Cross sections at 
energies away from the Z resonance peak have a much stronger dependence 
on higher order photonic corrections and their inclusion is of numerical im- 
portance at LEP2. Several improvements in this respect have been performed 
since ZFITTER version 6.21 and will be described in the following subsections. 
They have also been discussed in [6]. 



2. 1 Second order initial state fermion pair corrections 



Since ZFITTER version 6.21 there is an improved treatment of second order 
corrections for angular distributions and ApB- These corrections were applied 
in the leading logarithmic approximation as described in Ref. [2]. The option 
is accessed by a new flag: 

FBHO = - old treatment with photonic corrections only, 
FBHO = 1 - leading log fermion pair corrections are added. 
The non-singlet and singlet fermion pair contributions to the electron (positron) 
structure function of [8] (see Eq. (11) in [9] and Eq. (47) in [8]) are used (with 
different options according to the values of the flags ISPP, IPFC and IPSC). 
They are directly added to the corresponding photonic contributions [10, 11], 
which are governed by the flag FUNA. For pair corrections in general, the new 
value of flag IPTO = — 1 is added. It allows to calculate pure virtual pair con- 
tributions separately. This option can be used for comparisons or in the case, 
when the contribution of real pair emission is taken from another program. 
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2.2 Second order initial state QED corrections with acceptance cuts 



Since ZFITTER version 6.30, the second order initial state radiation (ISR) QED 
corrections in presence of angular cuts are improved [12]. A new option gov- 
erned by a new flag FUNA is implemented, with: 
FUNA=0 - old treatment, 
FUNA=1 - new treatment. 

The corrected treatment of these corrections takes the angular acceptance cuts 
ANGO, ANGl into account. The corrections are relevant for the angular distribu- 
tion and for the integrated forward-backward asymmetry. We use the leading 
logarithmic approximation by means of the electron structure function formal- 
ism [8,9]. In fact the differential angular distribution of the electron-positron 
annihilation process can be represented in a form analogous to that of the 
Drell-Yan process: 



= /' dx,V{xu Le) C dx2V{x2, L,)^^M je(5 _ (2.1) 
dc JO Jo dc 



with c = cos and 



s = XiX2S, (2.2) 

4:XlX2 

[Xi +X2- c{Xi - X2)Y 



J= ,^ , . 12 (2.3) 



where the structure functions T>{xi, L^) give the probability to flnd an electron 
(positron) with a reduced energy fraction xi(2) in the initial electron (positron). 
The Born-level annihilation cross section da{s, c)/dc is deflned in the center- 
of-mass reference frame of electron and positron with reduced energy fractions. 
The Jacobian J' is coming from the relation to the angles in the laboratory 
reference frame. The structure functions are taken in the leading logarithmic 
approximation keeping the first and second order photonic contributions: 



Vix, Le) = 5(1 - x) + ^(L, - l)pW(x) + I (|^(L, - 1))' P(2)(x),(2.4) 



with 
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mt 



lim{5(l-.)(: 



21nA + 



2 In A + - 



6(1 - A -x 



27r2 



l + x' \ 
l-x ' 



(2.5) 
(2.6) 

(2.7) 



+ 6(1- A-x)2 
1 + X 



( 2 ln(l - x) - In X + - 



+ 



Inx — 1 + X 



In order to avoid a double counting, we expanded formula (2.1) in a and take 
only the 0{a^) terms, which are then added to the full first order corrections. 
Where possible, we performed integrations over the angle and the energy frac- 
tions analytically. This allows to get a relatively fast code for these corrections. 
The calculation is realized with subroutine funang.f. It is compatible with 
the use of ICUT = 1,2,3. We just mention that for the contribution to the total 
cross section the complete 0{a'^) formulae are used. 



2. 3 Higher order photonic corrections from the initial- final state interference 



The exponentiation of photonic initial-final-state interference corrections was 
implemented according to Ref. [13]. It allows to take into account the most 
significant part of higher order corrections coming from the initial-final state 
interference. The corrections are relevant for the angular distribution and for 
the integrated forward-backward asymmetry. A combination of the one-loop 
initial-final interference corrections with the corresponding higher order effects 
from the exponentiation can be computed in ZFITTER by using flag INTF = 2. 
The old options are: 

INTF = 0: the initial-final state interference in photonic corrections is omitted; 
INTF = 1: the interference is taken in the one-loop approximation. 
The numerical effect of the exponentiation was discussed in Ref. [6]. It was 
found to be close to the one obtained by a slightly different exponentiation 
procedure [14]. 



2.4 Final state fermion pair production corrections 



Since ZFITTER version 6.30, the corrections from final state radiation (FSR) 
of fermion pairs are implemented [12], according to the formulae given in [15]. 
The option is governed by a new flag: 
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FSPP = - without FSR pairs, 
FSPP = 1 - with FSR pairs, additive, 
FSPP = 2 - with FSR pairs, muhiphcative. 

For the best approximation, the FSR pair corrections should be treated multi- 
phcative (FSPP = 2) with respect to the ISR photonic corrections (see Eq. (2.1) 
in [12]). The additive treatment of FSR pairs (FSPP = 1) can be used for a 
comparison. For the FSPP corrections, a cut on the invariant mass of the sec- 
ondary pair is accessible. In order to accommodate this cut value, the variable 
SIPP of the 

SUBROUTINE ZUCUTS ( INDF , ICUT , ACOL , EMIN , S_PR , ANGO , ANGl , SIPP) 

is now used. The dependence on this cut for realistic event selections was 
shown to be rather weak [12]. Therefore, the meaning of variable SIPP has 
been changed. Now it has nothing to do with cutting of the initial state pairs. 
In fact, there is no possibility to directly cut the initial state secondary pairs. 
The primary pair invariant mass cut (S_PR) is taken into account in the phase 
space of secondary pairs. 



3 Improved Born approximation for e+e" 
6.34 



vv in ZFITTER version 



For a study of the reaction 



e e 



one needs the effective Born approximation for 



(3.1) 



da 



dcosi!) 



i=e,fj,,T 



da{e^e — > ViVi 
dcos'd 



3 (7^ + ast + (yf 



(3.2) 



The specific property of this reaction is due to the interference of s channel Z 
boson and t channel W boson exchange for Vf.Ve production. For an application 
in KKMC [16], the corresponding formulae have been derived [17] from a 
related study [18,19] and implemented in ZFITTER version 6.34 (05 Feb 2001). 
The changes in the program consist of a modification of the weak charged 
current form factor (variable ROW) in subroutine coscut in zfbib6_34.f . The 
addition is ^ : 



a 



Scc^m = QED_CC-QED_NC = ^Ql 



^In^ 

2 s 



+ - In - 

2 s 



4 Li2(l) + 2 



(3.3) 



^ In the Fortran code until version 6.42, the factor of 4 at the r.h.s. of Equation 
(3.3) is lacking. 
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The variables SIGST and SIGT are influenced by this, they correspond to ast 
and CTj. Depending on the setting of flag ENUE there, part or all of them con- 
tribute to the prediction of da/dcos'd; see also Section B.2. 

In ZFITTER version 6.42, the file zfbib6_34.f is replaced by zfbib6_40.f. 
The latter file contains a subroutine cos cut where the described features are 
not accessible. The file zf EENN_34 . f is no longer part of the distribution. 



4 Atomic parity violation 



The global precision tests of the Standard Model may include an experimen- 
tal input from atomic parity violation measurements in heavy atoms. The 
observable quantity of interest is the so-called weak charge Qw'- 



Qw{Z, A) = -2 [{2Z + N) Ciu + {Z + 2N) Cm] • (4.1) 

Here, Z and are the numbers of protons and neutrons in the nucleus. 
The weak couplings involved are those which parameterize the electron-quark 
parity-violating Hamiltonian at zero momentum transfer: 



Hpv = -^{Ciue'j^'j^eu'y^u + C2ue'j^,eu-f^'j5U + CidC'jf.'j^ed'jf.d 

+ C'2de7^eJ7^75d). (4.2) 

In Born approximation, Ci„ = — ^(1—8/3 sin^ 6*1^/) and Ci^ = |(1— 4/3 sin^ 6'^/); 
generally, Ciq = 2aeVq and C2q = 2veaq, q = u,d. In [20], using the results 
of [21], the higher order corrections to atomic parity violation in the on-mass- 
shell renormalization scheme have been derived and the corresponding expres- 
sions were used for ZFITTER version 6.34 (26 Jan 2001) onwards. The formulae 
given in [22] were also reproduced. 

In subroutine ZU_APV the corresponding expressions for C\u^ C*ic/) C*2ji) are 
calculated. In fact, the sign conventions have been changed compared to [20] 
in order to share the definitions of the particle data group ^ . A sample use is 
prepared with subroutine ZF_APV, with N = 78 and Z = 55 for Caesium. Flag 
TUPV may be used for a study of the theoretical uncertainty, the default value 
is TUPV=1. 



^ We just mention that the signs of Equations (4) in [20] are not in accordance with 
Equation (3) there, while those of Equation (14) are. 
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5 Higher order electroweak corrections to Standard Model observ- 
ables 



In ZFITTER version 6.42, all known two- and three loop corrections of the 
Standard Model observables M^, sin^6'gg and Vz have consistently been in- 
cluded. The corresponding changes began with version 6.33, with subsequent 
additions until version 6.42. They are described to some detail in this section. 



5.1 The W boson mass 



ZFITTER, version 6.36 (21 July 2001), contained only the leading and next- 
to- leading corrections to Ar, obtained through an expansion in the heavy top 
quark mass. They were applied according to [23-28] and have been imple- 
mented in the package m2tcor . f . The W boson mass was then evaluated with 
the relation: 



1 + 



Ana 



V2MlG^il 



Ar) 



(5.1) 



by an iterative procedure, since Ar also depends on M^^. More accurate results 
on the electroweak corrections to muon decay, which successively appeared 
in [29-38], were found in terms of exact one- dimensional integral representa- 
tions at the order 0{a^). However, since the computation of general two loop 
integrals is rather slow, only fitting formulas as published in [32,37] were im- 
plemented in ZFITTER. Furthermore, since typical resummation prescriptions 
for Ar (see [39] and references therein) are problematic once the complete 
two-loop contributions are included, the new result has been directly applied 

to My.. 



With the option AMT4 = 5 the complete fermionic two loop corrections to the 
W boson mass [32] are used. Here the fermionic corrections encompass all 
two-loop contributions with at least one closed fermion loop. In addition to 
electroweak corrections, this setting also includes the leading 0{aag) [40-46] 
and next-to-leading ^(aa^) [47-49] QCD corrections. Contrary to the settings 
AMT4 < 4, the (9(aa^) corrections are incorporated exactly according to [49]. 
Previously, only the leading term in mf was used. 

The assignment AMT4 = 6, which is used from Summer 2004 onwards, enables 
the calculation of including complete (fermionic and bosonic) 0{a^) cor- 
rections [34,37], Oiaa^) [40-46] and Oiaa]) [49,50] QCD corrections, as well 
as leading three-loop corrections in an expansion in mf of order 0{a^) and 
0{a'^ag) [51]. These contributions were implemented into a fitting formula of 
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the following form [37]: 



M^ = Ml- ci dH - C2 dH^ + cg dH^ + C4(dh - 1) - cg da + ce dt (5.2) 
— C7 dt^ — cs dH dt + Cg dh dt — cio das + cn dZ, 



where 



dH 

dt = 
da : 



log { ) 
^ VlOO GeVy 



174.3 GeV 
Aa 



0.05907 



- 1, 



dh 



dZ 



da. 



VlOO GeV 
91.1875 GeV 



0.119 



- 1, 



(5.3) 



and the coefficients M^, ci, . . . , cn take the following values 



Ml 



C3 
Cg 



80.3799 GeV, ci 

0.0000890 GeV, C4 

0.5256 GeV, c^ 

0.0000659 GeV, cio 



0.05429 GeV, ca 

0.000161 GeV, C5 

0.0678 GeV, Cg 

0.0737 GeV, cn 



0.008939 GeV, 
1.070 GeV, 
0.00179 GeV, 
114.9 GeV. 



(5.4) 

With eq. (5.2) the full result for Myy is approximated to better than 0.5 MeV 
over the range of 10 GeV < < 1 TeV if all other experimental input values 
vary within their combined 2a region around the central values, as being used 
in eq. (5.3). 



For the new options AMT4 = 5,6 the algorithm used by ZFITTER for the esti- 
mation of the theory error due to the calculation of M^, based on applying 
different resummation procedures (see eq. (2.104) of [2]) cannot be used any- 
more, since the next term in the rrit expansion, which is of order OlGfj^m^M"^), 
is now known exactly. In contrast, for these options the uncertainties from 
unknown higher order contributions in the calculation of have been esti- 
mated in [32,37] and are set here to ±5 MeV and ±4 MeV, respectively. They 
can be simulated by varying the flag DMWW between —1 and 1, which is the 
only relevant flag to simulate theoretical uncertainties for for AMT4 = 5,6, 
whereas EXPR, IHIGS, HIG2, SCRE and SCAL are ignored. Furthermore, the use 
of AMT4 = 5, 6 is not compatible with the setting IMOMS > 1. 
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5.2 The effective leptonic weak mixing angle 



The amplitude for the decay of the Z boson into a pair of fermions is param- 
eterized in ZFITTER as: 



V,'" {Ml) = (z {2n)') I ^V2G,Mlsj7zlfl, [(1 + 75) - ^Qf\slA 

(5.5) 

This formula allows to define the effective weak mixing angle as: 

sin2^fg = Re(4)4, (5.6) 

where the coefficients p{ and k{ are called effective couplings of Z-decay. 

Complete fermionic two-loop corrections to the leptonic effective weak mixing 
angle sin^ 6^^^^ [52-54] are included with the option AMT4 = 6. They improve 
the realizations for AMT4 < 5 by including all electroweak two-loop contribu- 
tions with at least one closed fermion loop and go beyond the leading [55-58] 
and next-to- leading [59] terms in an expansion in m^. The corrections have 
been implemented directly using the numerical fitting formula as published 
in [52] , thereby assuring a fast evaluation: 

sin^ e'^l^ = So + d^Lu + d2L]j + d^L% + d^{^l - 1) + 4 A„ 

+4At + d^^l + 4At(A^ - 1) + rfgA,^ + dioA^, (5.7) 

with 



100 



^^-^^HWg^J' ^^-Toog^' 

Aa = 777777^ - 1' Ai= — — -4— - -1, (5. 



and 



0.05907 V178.0 GeV 

0.117 ' ^ 91.1876 GeV 



So = 0.2312527, di = 4.729 x 10~\ d2 = 2.07 x 10-^ 

4 = 3.85 X 10-^ c/4 = -1.85 X 10-^ d5 = 0.0207, 

4 = -0.002851, dr = 1.82 x 10~^ ds = -9.74 x 10-^ 

dg = 3.98 X 10-^ c/10 = -0.655. 



(5.9) 



In addition to the electroweak one- and two-loop corrections, this formula 
also includes two-loop C(aa<,) [40-46] and three-loop 0{aa'l) [49,50] QCD 
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corrections, as well as leading three-loop corrections in an expansion in 
of order 0{a^) and 0{a'^ag) [51]. In contrast to the settings AMT4 < 6, the 
three-loop QCD corrections are incorporated exactly according to [49], which 
goes beyond the leading term in used previously. Equation (5.7) reproduces 
the exact calculation with maximal and average deviations of 4.5 x 10~^ and 

1.2 X 10~^, respectively, as long as the input parameters stay within their 2a 
ranges and the Higgs boson mass in the range 10 GeV < < 1 TeV. 

With the new option AMT4 = 6 the theory error in the calculation of sin^ 9^^^^ 
is no longer obtained through different resummation formulas (see eq. (2.134) 
of [2]). Instead, the value of ±4.9 x 10~^, as estimated in [52], is used. For 
this case the theory uncertainty can be simulated by varying the flag DSWW 
between —1 and 1. 

5.3 Partial and total Z boson widths 

The Z boson decay width is defined through the effective couplings, /t{ and 
p{, of eq. (5.5). The form factor is obtained from the effective weak mixing 
angle sin^ ^eff**, as presented in subsection 5.2. For Z boson decays into a pair 
of fermions f,fj^l, the two-loop corrections to k{ are currently implemented 
in an approximate way only, while for T{Z ^ bb), no two- loop corrections 
beyond the leading mf term are available. For the form factor p{, the two- 
loop corrections beyond the next-to-leading mf expansion are still missing. 
This quantity is therefore computed identically for the choices AMT4 = 4, 5, 6. 

5.4 Electroweak form factors for e+e^ ff 

The radiative corrections to the cross-sections and asymmetries for the pro- 
cess e~^e~ — > // are parameterized in ZFITTER by electroweak form factors, 
Pef, He) i^f and Ke/, defined through a Z-boson exchange amplitude: 




-MQe\slrKe{s,t)-ff, (g) 7^,(1+75) - 4\Q f\sl,Kf{s,t)-f^{l + 75) 7^ 




(5.10) 



where 



Xz{s) 



G, Ml s 
1/2 Sttq; s — m| 



and 




iM^T^is). 



(5.11) 
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In the leading pole approximation, the following relations hold: 

Pef = \J P'z Pz^ ~ '^Z5 ~ ^ef ~ ^z^z' (5.12) 

However, generally, pe/, ^^e, t^f, and Kg/ also include ■y-Z interference effects, 
corrections from non-resonant Z and 7 exchanges, and non-factorizable box 
contributions. With the option AMT4 = 6 the coefficients Ke and k/, for / 7^ 6, 
are calculated with all known two- and three-loop effects, as discussed in the 
previous sections. The corresponding factorizable vertex corrections to Kg/ 
are incorporated similarly. The electroweak box contributions to the form 
factors include the complete one-loop order correction, which is sufficient for 
present precision. For the form factor p^f, complete two- loop corrections are 
still missing, and it is computed including two-loop corrections of the order 
0{GlmlMl) only for all options AMT4 = 4, 5, 6. 

For the Z ^ bb channel, no two-loop electroweak corrections beyond the lead- 
ing mf terms are available so far. Furthermore, in ZFITTER versions before 
6.42, a mismatch occurred in the treatment of the bb final state. The subrou- 
tine GDEGNL, which computes the leading two-loop corrections, is called at two 
places in the program: (z) in subroutine ZWRATE the effective couplings intro- 
duced in section 5.2 are used for the computation of Z partial widths; (ii) 
the interfaces ZUTHSM, ZUTPSM, ZULRSM and ZUATSM, on the other hand, calcu- 
late the cross-sections and asymmetries from the weak form factors pef{s,t), 
Ke{s,t), Kf{s,t) and Kef{s,t) in the subroutine ROKANC. For INDF = 9, i.e. the 
bb final state, all these form factors were calculated in one-loop approximation 
with the leading two-loop term. As a result, in ROKANC the initial-state 
Ze~^e~ form factors for all other final states are generated including next- 
to-leading two-loop corrections, while for the bb final state these corrections 
were missing. This mismatch also affects the ZFITTER interfaces ZUXSA, ZUTAU 
and ZUXSA2, which use the language of effective couplings [2], since they are 
defined to coincide as close as possible with the complete Standard Model 
prediction in ROKANC if the effective couplings coincide with their Standard 
Model analogue. 

The problem has been alleviated since ZFITTER version 6.41 (15 October 2004). 
In contrast to the older implementations, K,e{s,t) and Kf{s,t) are not treated 
symmetrically anymore for INDF = 9, but two-loop electroweak corrections 
are included in Ke{s,t) for AMT4 > 4, not yet in Kb(s,t). The treatment of 
Pef{s,t) and Kef{s,t) has been changed accordingly^. Here one can use the 
fact that the presently known two-loop contributions factorize into initial-state 
and final-state corrections. For a more detailed discussion see [60]. 



^ In ZFITTER version 6.41 an error occured in the treatment of a form factor 
Pef{s,t). It affected the bb cross-section only, and was corrected in ZFITTER ver- 
sion 6.42. 
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6 Input parameters and pseudo observables 



For the calculations of pseudo observables (POs), such as partial Z decay 
widths or effective coupling constants, ZFITTER uses the DIZET package which 
calculates such quantities employing the on-mass- shell ( OMS) renormalization 
scheme [61,62] within the Standard Model. The user has to provide a set of 
values for the so-called input parameter set, which are then used by ZFITTER 
and DIZET to calculate the POs. A standard set of imput parameters is given 



• The electromagnetic coupling constant at the Z pole, or, more precisely, its 
shift due to the 5-quark flavour hadronic vacuum polarisation, Aa|^^|j at the 
Z pole; 

• The strong coupling constant at the Z pole, as] 

• The pole masses of Z boson, top quark and Higgs boson. 

These variables are the main physics values to be provided to ZFITTER and 
DIZET by the user. The Fermi constant is treated as a constant in the pro- 
gram (see flag GFER). In order to use the user-supplied value of the hadronic 
vacuum polarisation, flag ALEM=2 must be used. In order to use the latest set 
of electroweak radiative corrections (see [37,52] and Section 5), flag AMT4=6 
must be used. 

Selected POs calculated by ZFITTER and DIZET are listed in Tab. 6.1, namely: 
(i) the W boson mass and the on-shell electroweak mixing angle = 1 — 
M^/M"^; (ii) the partial Z decay widths, including the invisible width, simply 
equal to SFj^, and the the total hadronic width, equal to the sum of the five- 
flavour quarkonic widths, and the total width; (iii) ratios of widths and pole 
cross-sections as used by, e.g., the LEP EWWG; (iv) the effective electroweak 
mixing angle sin^ 6'fg according to Eq. (5.6) and the pf parameter pf = Re{pz), 
for leptons and heavy fermions b and c; and (v) the asymmetry parameter 
Af and the forward-backward pole asymmetries A^-^ for leptons and heavy 
fermions b and c. 

The partial decay widths are defined inclusively, i.e., they contain all real and 
virtual corrections. The ratios for leptons / and quarks q are defined as: 



by: 





(6.1) 



while the hadronic and leptonic pole cross-sections are defined as: 



127r 



TeTh 

M|F| 



.0 



127r 



M|F| ■ 



(6.2) 
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The complex variable (7^ is defined as the ratio: 

gi = ^ = l-4\Qf\4sl (6.3) 

of the complex effective vector and axial couplings of the Z boson to the 
fermion /: 



Vf = ^piiTi -2Qf4sl), (6.4) 

«/ = ^/7zTi. (6.5) 
Hence one may write: 

Re{gf) = l-A\Qf\sm'9ls. (6.6) 

The complex form factors p;^, k;^ are defined through the amplitude of the Z 
boson decay into a pair of fermions, as in eq. (5.5). 

With these definitions, the asymmetry parameters are then given in terms of 
the real part of g^'- 

Re 

Af = 2 (6.7) 
l+(Reglj 

The forward-backward pole asymmetries are mere combinations of these asym- 
metry parameters: 

Kl = l^eAj. (6.8) 

In terms of the real part of the complex p;^ parameter, pf = Re(p;^), and 
the (real) effective electroweak mixing angle sin^ O^^ defined earlier, the real 
effective vector and axial-vector coupling constants as quoted by the LEP 
EWWG are then defined as: 

g^f = (6.9) 
gvf = ^{Ti -2Q J sin' 9l^). (6.10) 

With these definitions, it follows that: 

^ = Mgl), (6.11) 

9Af 

so that the asymmetry parameters are then equivalently given as: 

Af = (6.12) 

l + igvf/gAfV 
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Observable 


Value 


M„ [GeV] 


80.3613 




0.22335 


[MeV] 


167.219 


Te [MeV] 


83.990 


[MeV] 


83.990 


[MeV] 


83.800 


Tc [MeV] 


299.966 


Tfe [MeV] 


375.729 


Tinv [GeV] 


0.501657 


T/, [GeV] 


1.741507 


[GeV] 


2.494944 



Table 6.1 



Observable 


Value 




20.73462 


Rl 


0.215750 


R'c 


0.172245 


al [nb] 


41.4826 


a? [nb] 


2.0006 


singer 


0.231548 


sin^ el^ 


0.233032 


sin2 


0.231442 


Pe 


1.005165 


Pb 


0.993943 


Pc 


1.005860 



Observable 


Value 




0.146813 


Ab 


0.934554 




0.667779 




0.016166 




0.102904 


^1'^ 


0.073529 



Predictions for pseudo observables calculated with ZFITTER and DIZET for Aa^J^ = 
0.02758, as = 0.118, Mz = 91.1875 GeV, Mt = 175 GeV and Mh = 150 GeV, 
and flags ALEM=2 and AMT4=6. 



6.1 Pseudo observables in common blocks o/ DIZET 



The channel dependent quantities are stored in common block COMMON/ CDZRKZ/ 
and in array PARTZ(0: 11): 



CDMMDN/CDZRKZ/ARRDFZ (0:10), ARKAFZ (0:10), ARVEFZ (0 : 10) , ARSEFZ (0:10) 
& ,ARDTFZ(0:10) ,AIR0FZ(O: 10) , AIKAFZ (0 : 10) ,AIVEFZ(0 : 10) 

* 

DIMENSION NPAR(25) ,ZPAR(30) ,PARTZ(0:11) ,PARTW(3) 
* 

SUBROUTINE DIZET (NPAR,AMW ! NPAR : FLAGS; AMW : INPUT/OUTPUT 
& ,AMZ,AMT,AMH,DAL5H,V_TBA,ALSTR ! INPUT 

& ,ALQED,ALSTRT,ZPAR,PARTZ,PARTW) ! OUTPUT 



The correspondences have not been changed compared to Section 2.5.1 of [2]: 
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ARR0FZ(0:10) 
AR0TFZ(0:10) 
ARKAFZ(0:10) 
ARVEFZ(0:10) 
AIR0FZ(0:10) 
AIKAFZ(0:10) 
AIVEFZ(0:10) 
ARSEFZ(0:10) 



{pi)' , 



6.13) 
6.14) 
6.15) 
6.16) 
6.17) 
6.18) 
6.19) 
6.20) 
6.21) 



Re pi , 

Imp{, 
Im k{ , 




PARTZ(0:11) 



The usual ZFITTER channel assignments as given in Figure 5 of [2] are used, and 
all quantities but (p{) have been introduced aheady. The (p{) is discussed in 
Section 2.5.1 of [2]. We mention also here that both options MISC=1,0 (with 
(p{) or with |p{|) are used in the Model Independent interfaces of ZFITTER, 
see Appendix C. 

Note that since the convention of the LEP EWWG is to use the real parameter 
p/ = Re(p^), not {pz)\ the array AROTFZ is generally used, while ARROFZ, 
RENFAC and SRENFC are usually ignored. 



7 Subroutine ZFTEST 

The ZFITTER distribution includes subroutine ZFTEST. With ZFTEST the user 
may test whether ZFITTER has been properly installed. The subroutine cal- 
culates cross-sections and asymmetries as functions of ^/s near the Z peak, 
below, and above. A sample file zfmai6_42.f runs ZFTEST: 



CALL ZFTEST (0) 
END 

The numerical output should reproduce the Tables given in Subsection 7.1. 
7.1 ZFTEST results 

This Section contains the standard test-outputs, produced by a call to ZFTEST (0) . 
The argument of ZFTEST sets the flag MISC. The default value is MISC=0. The 



* 



* MAIN to call ZFTEST with version 6.42 



* 
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user's result should be identical to the sample output, apart from a possible 
flip in the last digit. 
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** This is ZFITTER version 6.42 

** 05/05/18 

** http : //www . ifh. de/theory/publist . html 
ZUINIT> ZFITTER defaults: 



ZFITTER 


flag 


values : 












AFBC 


1 


SCAL 





SCRE 





AMT4 


4 


BORN 





BDXD 


1 


CONY 


1 


FIHR 


1 


FDT2 


3 


GAMS 


1 


DIAG 


1 


INTF 


1 


BARB 


2 


PART 





PDWR 


1 


PENT 





ALEM 


3 


qCDC 


3 


VPDL 


1 


WEAK 


1 


FTJR 


1 


EXPR 





EXPF 





HIGS 





AFMT 


3 


CZAK 


1 


PREC 


10 


HIG2 





ALE2 


3 


GFER 


2 


ISPP 


2 


FSRS 


1 


HISC 





MISD 


1 


IPFC 


5 


ipse 





IPTO 


3 


FBHD 





FSPP 





FDNA 





ASCR 


1 


SFSR 


1 


EHUE 


1 


TUPV 


1 


DMWW 





DSUW 






















ZFITTER cut values : 



to 



IKDF 


ICUT 


ACOL 


EMIN 


s 


_PR 


ANGO 


ANGl 


SPP 










00 





00 





00 





00 


180 


00 


0.00 


1 







00 





00 





00 





00 


180 


00 


0.00 


2 







00 





00 





00 





00 


180 


00 


0.00 


3 







00 





00 





00 





00 


180 


00 


0.00 


4 







00 





00 





00 





00 


180 


00 


0.00 


B 







00 





00 





00 





00 


180 


00 


0.00 


6 







00 





00 





00 





00 


180 


00 


0.00 


7 







00 





00 





00 





00 


180 


00 


0.00 


8 







00 





00 





00 





00 


180 


00 


0.00 


9 







00 





00 





00 





00 


180 


00 


0.00 


10 







00 





00 





00 





00 


180 


00 


0.00 


11 







00 





00 





00 





00 


180 


00 


0.00 



ZFITTER input parameters: 
DAL5H - 0.0280398093 
ALqED5= 128.886183 
ZMASS - 91.18760; THA 
HMASS - 100.00000 
DAL5H - 0.02804; 
ALFAS = 0.11700; 



ALqED5 - 128.88618 
ALFAT - 0.10637 



ZFITTER intermediate results: 
WMASS - 80.39522; SIN2TU - 0.22270 

ALPHST - 0.11700; 

qCDCOR - 1.00000 1.03938 1.04639 1.03877 1.03177 1.03941 1.04611 1.03877 1.03177 1.03938 1.046 



nu.nubar 


167 


299 


1 


008546 


1 


008546 





231036 


e+,e- 


84 


037 


1 


005790 


1 


005643 





231417 


mu+.mu- 


84 


037 


1 


005790 


1 


005643 





231417 


tau+,tau- 


83 


847 


1 


005790 


1 


005643 





231417 


u.ubar 


300 


128 


1 


006390 


1 


006339 





231310 


d.dbar 


383 


033 


1 


007315 


1 


007306 





231183 


c,cbar 


300 


059 


1 


006390 


1 


006339 





231310 


s,sbar 


383 


033 


1 


007315 


1 


007306 





231183 


t,tbar 





000 





000000 





000000 





000000 


b,bbar 


375 


593 





993885 





993885 





233004 


hadron 


1741 


846 














total 


2495 


664 















1.03937 1.02506-0.00002 0.19863 



SQRT(S) = 89.1875992 



< Cross Section > < Asymmetry > <~Tau__Pol~> < A_LR > 



INDF 


ZUTHSH 


ZUXSEC 




ZDXSA 


ZUXSA2 


ZUXAFB 


ZDTHSM 


ZUXSA 


ZUXSA2 


ZUXAFB 


ZUTPSH ZDTAD 


ZDLESM 


ZUALR 








.79464 





.79464 
































1 





.41291 





.41291 





.41291 


0.41291 


0.41291 


-0. 


1886- 


-0.1886- 


-0. 


1885- 


-0. 1885 












2 





.41309 





.41309 





.41309 


0.41309 


0.41309 


-0. 


1886- 


-0.1886- 


-0. 


1885- 


-0. 1885 












3 





.41223 





.41223 





.41223 


0.41223 


0.41223 


-0. 


1888- 


-0.1888- 


-0. 


1888- 


-0. 1888- 


-0. 1206-0. 1206 










4 


1 


.44131 


1 


.44131 


1 


.44131 






-0. 


0473- 


-0.0473 













.0841 





.0000 


5 


1 


.81653 


1 


.81653 


1 


.81663 






0, 


0544 


. 0544 













.1031 


2 


. 1932 


6 


1 


.44105 


1 


.44105 


1 


.44105 






-0. 


0474- 


-0.0474 













.0841 


2 


.1932 


7 


1 


.81654 


1 


.81654 


1 


.81663 






0, 


0544 


. 0644 













.1031 


2 


.0196 


8 





. 00000 





. 00000 





, 00000 






0, 


0000 


. 0000 













.0000 


2 


. 1951 


9 


1 


.78022 


1 


.78022 


1 


, 78022 






0, 


0542 


. 0642 













.1026 





.0000 


10 


8 


.29565 


8 


.29565 



























.0963 





.0000 


11 





.70082 





.70101 





. 70090 


. 70090 




0, 


4313 


0.4312 





.4312 















SQRTCS) = 91.1875992 



< Cross Section > < Asymmetry > < — Tau_Pol — > < A_LR > 



INDF 


ZUTHSH 


ZUXSEC 




ZUXSA 


ZUXSA2 


ZUXAFB 


ZDTHSM 


ZUXSA 


ZUXSA2 


ZUXAFB 


ZUTPSM ZDTAD 


ZDLRSM 


ZUALR 





2 


.91319 


2 


.91319 
































1 


1 


.47692 


1 


.47692 


1 


. 47692 


1 


. 47692 


1 . 47692 





.0003 





.0003 


0, 


0003 


. 0003 








2 


1 


.47745 


1 


.47745 


1 


. 47746 


1 


. 47746 


1 . 47746 





.0003 





.0003 


0. 


0003 


. 0003 








3 


1 


. 47438 


1 


. 47438 


1 


. 47438 


1 


. 47438 


1 . 47438 





.0001 





.0001 


0, 


0001 


0.0001- 


-0.1431-0.1431 






4 


5 


.24201 


5 


.24201 


6 


.24201 











.0615 





.0616 










0.1418 


.0000 


6 


6 


.67109 


6 


.67109 


6 


.67109 











.0971 





.0971 










0.1441 2 


.1932 


6 


5 


.24108 


5 


.24108 


6 


.24109 











.0616 





.0616 










0.1418 2 


. 1932 


7 


6 


.67112 


6 


.67112 


6 


.67113 











.0971 





.0971 










0.1441 2 


.0196 


8 





. 00000 





.00000 





. 00000 











.0000 





.0000 










0.0000 2 


. 1951 


9 


6 


.54153 


6 


.54153 


6 


.64163 











.0982 





.0982 










0.1444 


.0000 


10 


30 


. 36683 


30 


.36683 




























0.1434 


.0000 


11 


1 


.37174 


1 


.37182 


1 


.37174 


1 


.37174 







. 1869 





.1869 


0. 


1869 











SQRT(S) = 93.1875992 



< Cross Section > < Asymmetry > < — Tau_Pol — > < A_LR > 



INDF 


ZUTHSH 


ZUXSEC 




ZDXSA 


ZUXSA2 


ZUXAFB 


ZDTHSM 


ZDXSA 


ZUXSA2 


ZUXAFB 


ZUTPSM ZDTAD 


ZDLRSM 


ZUALR 





1 


.22102 


1 


.22102 




























1 





.62880 





.62880 





. 62880 


. 62880 


. 62880 


0, 


1209 





.1209 


0. 1209 


0.1209 








2 





.62905 





.62905 





. 62906 


. 62906 


. 62906 


0, 


1209 





.1209 


0. 1209 


0.1209 








3 





.62781 





.62781 





.62781 


0.62781 


0.62781 


0, 


1208 





.1208 


0. 1208 


0.1208- 


-0.1630-0.1630 






4 


2 


.21270 


2 


.21270 


2 


.21270 






0, 


1293 





,1293 








0.1763 


.0000 


6 


2 


.80592 


2 


.80592 


2 


.80692 






0, 


1235 





.1236 








0.1696 2 


. 1932 


6 


2 


.21234 


2 


.21234 


2 


,21236 






0, 


1295 





,1296 








0.1763 2 


.1932 


7 


2 


.80593 


2 


.80593 


2 


.80694 






0, 


1235 





,1235 








0.1696 2 


.0196 


8 





.00000 





.00000 





. 00000 






0, 


0000 





,0000 








0.0000 2 


. 1951 


9 


2 


.75285 


2 


.75285 


2 


. 76286 






0, 


1253 





,1253 








0.1704 


.0000 


10 


12 


.78974 


12 


.78953 
























0.1720 


.0000 


11 





. 52345 





.52335 





. 62336 


. 62336 




0, 


2174 





,2174 


0.2174 











SQRTCS) = 100. 



< Cross Section > < Asymmetry > < — Tau_Pol — > < A_LR > 



INDF 


ZUTHSH 


ZUXSEC 




ZDXSA 


ZUXSA2 


ZUXAFB 


ZDTHSM 


ZUXSA 


ZUXSA2 


ZUXAFB 


ZUTPSM ZDTAU 


ZULRSM 


ZUALR 








.20718 





.20718 


































1 





.11707 





. 11707 


0. 


11707 





.11707 


0.11707 


0. 


2404 





.2404 





.2404 


. 2404 










2 





.11714 





. 11714 


0. 


11715 





.11716 


0.11714 


0. 


2403 





.2403 





.2403 


. 2403 










3 





. 11696 





. 11696 


0. 


11696 





.11696 


0.11696 


0. 


2405 





.2405 





.2405 


. 2405- 


-0.1520-0.1520 








4 





.39006 





.39006 


0. 


39006 








0. 


2003 





.2003 










. 2078 





.0000 


5 





.48082 





.48082 


0. 


48082 








0. 


1518 





.1518 










0.1962 


2 


.1932 


6 





.39002 





.39002 


0. 


39002 








0. 


2006 





.2006 










. 2078 


2 


. 1932 


7 





.48083 





.48083 


0. 


48083 








0. 


1518 





.1518 










0.1962 


2 


.0196 


8 





. 00000 





. 00000 
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8 Summary 



The description of the ZFITTER package has been updated to version 6.42. 
Besides a short introduction to the physics contents of the program additions, 
we reproduce, in Appendices, also the technical details like user flags and 
interface subroutines as complete as necessary for a convenient use of the 
program. 

ZFITTER covers most of the radiative corrections of a practical relevance in 
the foreseeable future. Yet, several corrections are still needed to close the two 
loop program of data analysis with ZFITTER: 

• cross section (or decay rate) asymmetries, with the Z boson coupling to 
light quarks or h quarks; 

• electroweak corrections to the Z boson width; 

• bosonic electroweak corrections to all asymmetries and the effective weak 
mixing angle. 

The inclusion of the two loop electroweak corrections to the Z\k) vertex will 
additionally require to account for the mass of the h quark in the one loop 
electroweak corrections. Nonetheless, it was checked that the one loop cor- 
rections with a massive h quark would give negligible effects for the accuracy 
reached at LEP. 

Some bug reports concerning ZFITTER versions between version 6.21 and ver- 
sion 6.42 will be found at the webpage [3]. 
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A DIZET user guide 

This Appendix describes technical details of the DIZET package. Not all of 
them have been influenced by the updates over the years. For the convenience 
of the user, we nevertheless decided to give a complete overview and will repeat 
a substantial part of the material, which was already presented in Section 4.1 
of [2]. 

A.l Structure o/ DIZET 

A first call of subroutine DIZET returns various pseudo-observables, the W- 
boson mass, weak mixing angles, the Z-boson width, the W^-boson width and 
other quantities. After the first call to DIZET, several subroutines of DIZET 
might be used for the calculation of form factors and couplings. This is de- 
scribed in Section 4.1.1 of [2]. 

A. 2 Input and output o/ DIZET 

The DIZET argument list contains Input, Output and Mixed (I/O) types of 
arguments: 



CALL DIZET (NPAR , AMW , AMZ , AMT , AMH , DAL5H , ALQED , ALSTR , ALSTRT , ZPAR , PARTZ , PARTW) 



A. 2.1 Input and I/O parameters to be set by the user 
Input: 

NPARCl : 25) , INTEGER*4 vector of fiags 
AMT = rrit - t-quark mass 
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AMH = Mh - Higgs boson mass 

ALSTR = as(M|) - strong coupling at s = M| 



I/O: 

AMW = Mw, W boson mass, input if NPAR(4) = 2,3, but is being calculated 
for NPAR(4) = 1 

AMZ = Z boson mass, input if NPAR(4) = 1,3, but is being calculated for 
NPAR(4) = 2 

DAL5H = Aa!f^(M|), hadronic vacuum polarization 

The M|, M^, and A«|f^(M|) cannot be assigned by a parameter statement 
(input / output variables) . 



A. 2. 2 Output of the TilZKY package 

ALQED = a(M|), calculated from Aajf^(Mj), see description of flag ALEM in 

Subsection A. 3 
ALSTRT = ag{mj) 

ZPAR(l) = DR=Ar, the loop corrections to the muon decay constant 
ZPAR(2) = DRREM = Arrem, the remainder contribution O {a) 
ZPAR(3) = SW2 = s^, squared of sine of the weak mixing angle defined by 
weak boson masses 

ZPAR(4) = GMUC = G^, muon decay constant, if NPAR(4) = 1,2, is set in 
CONSTl depending on flag NPAR(20), see Subsection A. 3. (It should be cal- 
culated if NPAR(4)=3 from Mz.M^, but then it will deviate from the ex- 
perimental value.) 

ZPAR(5-14) - stores effective sines for all partial Z-decay channels: 

5 - neutrino 

6 - electron 

7 - muon 

8 - r lepton 

9 - up quark 

10 - down quark 

11 - charm quark 

12 - strange quark 

13 - top quark (presently equal to up quark) 

14 - bottom quark 
ZPAR(15) = ALPHST = ag{Ml) 

ZPAR(16-30) = QCDC0R(0-14), QCDCOR(I) - array of QCD correction factors 
for quark production processes and/or Z boson partial width (channel i) 
into quarks. Enumeration as follows: 
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QCDCOR(O) = 1 
QCDCOR(l) = Ry{Ml) 
QCDC0R(2) = Rl{Ml) 
QCDCORO) = Ri{Ml) 
QCDC0R(4) = Ri{Ml) 
QCDC0R(5) = Ry{Ml) 
QCDC0R(6) = R'^j.^Ml) 
QCDC0R(7) = Ry{Ml) 
QCDC0R(8) = R\{Ml) 
QCDCORO) = R^M^) 
QCDCOR(IO) = i?"(M|) 
QCDCOR(ll) = Ri{Ml) 
QCDC0R(12) = R\{Ml) 
QCDCORdS) = -R:^-singlet vector correction 
QCDC0R(14) = /i, corrections to ApB [63] 



(A.l) 



foreseen for 
foreseen for 



tt- channel 
tt-channel 



PARTZ(I) - array of partial decay widths of the Z-boson: 
1=0 neutrino 
1=1 electron 
1=2 muon 
1=3 tau 
I = 4 up 
1=5 down 
1=6 charm 
1=7 strange 

1=8 top (foreseen, not realized) 
1=9 bottom 
I = 10 hadrons 
I = 11 total 

PARTW(I) - array of partial decay widths of the VT-boson^ for the channels: 
1 = 1 one leptonic 
1 = 2 one quarkonic 
1 = 3 totoal 



A. 3 The flags used by DIZET 



Since the DIZET package may be used as stand-alone in order to compute 
POs we present here a short description of all flags in DIZET. The flag values 

^ The calculation of the W width [64] follows the same principles as that of the Z 
width and is realized in subroutine ZWRATE of DIZET. Since the W width is not that 
important for the description of fermion pair production, we do not go into details. 
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must be filled in vector NPAR(1 : 25) by the the user. Most of these flags overlap 
with the flags set in user subroutine ZUFLAG called by ZFITTER, however, in the 
stand-alone mode ZUFLAG need not be called. We will show the correspondence 
between the flag names CHFLAG and the flag values IVALUE used inside DIZET, 
called with 

CALL ZUFLAG CHFLAG ' , IVALUE) . 

The description is given in the order of the vector NPAR(1:25). Flag values 
marked as presently not supported are not recommended. For instance, 
they may be chosen for backward compatibility with respect to earlier versions 
of the code. 

NPAR(l) = IHVP ZUFLAGCVPOL' ,IHVP) - Handling of hadronic vacuum 
polarization: 

IHVP = 1 (default) by the parameterization of [65] 

IHVP = 2 by effective quark masses of [66, 67] presently not 

IHVP = 3 by the parameterization of [68] supported 

NPAR(2) = IAMT4 ZUFLAG AMT4 ' , IAMT4) - Re-summation of the leading 
O {Gi^.m'l) electroweak corrections, see Section 5: 
IAMT4 = no re-summation 

IAMT4 = 1 with re-summation recipe of [69] presently 
IAMT4 = 2 with re-summation recipe of [43] not 
IAMT4 = 3 with re-summation recipe of [70] supported 
IAMT4 = 4 (default) with two-loop sub-leading corrections and re-summation 
recipe of [23-28] 

IAMT4 = 5 with fermionic two-loop corrections to according to [29,30, 
32] 

IAMT4 = 6 with complete two-loop corrections to My^,. [37] and fermionic 
two-loop corrections to sin^ 6*^^' [52] 

NPAR(3) = IQCD -> ZUFLAG QCDC ' ,IQCD) - Handling of internal QCD cor- 
rections of order O (aOg): 
IQCD = no internal QCD corrections 
IQCD = 1 by Taylor expansions (fast option) of [71] 
IQCD = 2 by exact formulae of [71] 
IQCD = 3 (default) by exact formulae of [42] 

NPAR(4) — IMOMS - Choice of two input/output parameters from the three 
parameters M^, M^^i/j: 

IMOMS = 1 (defauh) input G^.Mz] output M^y), see Eq. (5.1) 
IMOMS = 2 input G^,M^] output M^, 

IMOMS = 3 input M^^Myv] output G^, foreseen, not realized 

NPAR(5) = IMASS - Handling of hadronic vacuum polarization in Ar; for 
tests only: 
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IMASS = (default) uses a fit to data 
IMASS = 1 uses effective quark masses 

NPAR(6) = ISCRE ZUFLAG( 'SCRE' , ISCRE) - Choice of the scale of the 
two-loop remainder terms of Ar with the aid of a conversion factor /, for 
details see [2]. 

ISCRE = (default) scale of the remainder terms is -ft'scaie = 1 
ISCRE = 1 scale of the remainder terms is -ft'scaic = 

ISCRE = 2 scale of the remainder terms is -ft'scaic = 

J 

NPAR(7) = lALEM ZUFLAG ( ' ALEM ' , lALEM) - Controls the usage of a(M|), 
see flowchart in [2]. Inside DIZET, however, its meaning is limited: 
lALEM = or 2 Aa|f'*(M|) must be supplied by the user as input to the 
DIZET package 

lALEM = 1 or 3 Aa|f'*(M|) is calculated by the program using a parame- 
terization IHVP (default: IALEM=3) 

For details see the complete discussion about this flag in Sections 2.8 and 4.2.2 
of [2]. 

NPAR(8) = IMASK - Historical relict of earlier versions. Presently unused. 

NPAR(9) = ISCAL ^ ZUFLAG( 'SCAL' , ISCAL) - Choice of the scale of ag(^mj) 
ISCAL = (default) exact AFMT correction [47] 
ISCAL = 1,2,3 options used in [72], presently not supported 
ISCAL = 4 Sirlin's scale ^ = 0.248 [73] 

NPAR(IO) = IBARB ^ ZUFLAG( 'BARB IBARB) - Handling of leading C» (C^m 
corrections: 

IBARB = corrections are not included 

IBARB = 1 corrections are applied in the limiting case: Higgs mass negligi- 
ble with respect to the top mass, [74] 
IBARB = 2 (default) analytic results of [75] approximated by a polynomial 

[76] 

These options are inactive for AMT4 = 4. 

NPAR(ll) = IFTJR ZUFLAGCFTJR' ,IFTJR) - Treatment of O^G^a^ml) 
FTJR corrections [77], 
IFTJR = without FTJR corrections 
IFTJR = 1,2 with FTJR corrections (default 1FTJR=1) 

Inside DIZET its meaning is limited. See complete discussion about this flag 
in [2]. 

NPAR(12) = IFACR -> ZUFLAG( 'EXPR' , IFACR) - Realizes different expan- 
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sions of Ar [26,59]: 

IFACR = (default) realizes the so-called OMS-I renormalization scheme 
IFACR = 1 intermediate step from OMS-I to OMS-II renormalization scheme 
IFACR = 2 approaches the spirit of the OMS-II renormalization scheme, a 
fully expanded option 

NPAR(13) = IFACT ZUFLAG( 'EXPF' , IFACT) - To simulate theoretical un- 
certainties different expansions of the formfactors p and k are realized in 
complete analogy to the flag IFACR: 
IFACT = (default) OMS-I renormalization scheme 

IFACT = 1 intermediate step from OMS-I to OMS-II renormalization scheme 
IFACT = 2 approaches the OMS-II renormalization scheme 

NPAR(14) = IHIGS ^ ZUFLAGCHIGS' ,IHIGS) - Switch on/off resummation 
of the leading Higgs contribution: 

IHIGS = (default) leading Higgs contribution is not re-summed 
IHIGS = 1 leading Higgs contribution is re-summed 

NPAR(15) = lAFMT ^ ZUFLAG( ' AFMT' , lAFMT) - Includes the three-loop cor- 
rections AFMT~ 5°°' [47] (see also description of flag SCAL): 
lAFMT = without AFMT correction 
lAFMT = 1 correction 0{G frn^al) is included 

lAFMT = 2 corrections 0{Gfmlal) and O^CfM^al + \og{mf)) are in- 
cluded 

lAFMT = 3 (default) corrections 0{Gfmfa^^), 0{GfM^al + \og{m^)) and 
0{G fMz/mlal) are included 

NPAR(16) = lEWLC - Treatment of the remainder terms of p and k (used in 
ROKAPP together with obsolete option AMT4 = 1-3): 
lEWLC = all remainders are set equal to zero 
lEWLC = 1 (default) standard treatment 

NPAR(17) = ICZAK ^ ZUFLAG( 'CZAK' , ICZAK) - Treatment of the CKHSS non- 
factorizable 0{aag) corrections ^ew/qcd to the quarkonic width, Vq [78, 
79]: 

ICZAK = without CKHSS corrections 
ICZAK = 1 (default) with CKHSS corrections 

ICZAK = 2 with CKHSS for pseudoobservables and without CKHSS in case of 
realistic observables 

Inside the weak library DIZET the meaning of ICZAK is limited. See also dis- 
cussion about this flag in Subsection B.2. 

NPAR(18) = IHIG2 -> ZUFLAG( ' HIG2 ' , IHIG2) - Handling of the quadrati- 
cally enhanced two-loop Higgs contributions to Ar [80,81]: 
IHIG2 = without Higgs corrections 
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IHIG2 = 1 with Higgs corrections 



NPAR(19) = IALE2 ZUFLAG( ' ALE2 ' , IALE2) - Treatment of leptonic cor- 
rections to Aa: 

IALE2 = for backward compatibility with versions up to v.5.12 

IALE2 = 1 with one-loop corrections 

IALE2 = 2 with two-loop corrections [82] 

IALE2 = 3 (default) with three-loop corrections [83] 

NPAR(20) = IGFER ZUFLAG( 'GFER' , IGFER) - Handhng of QED correc- 
tions to the Fermi constant: 

IGFER = for backward compatibility with versions up to v.5.12 
IGFER = 1 one-loop QED corrections for Fermi constant [84-86] 
IGFER = 2 two-loop QED corrections for Fermi constant [87, 88] 

NPAR(21) = IDDZZ - Used in ZWRATE for internal tests: 
IDDZZ = RQCDV(A) are set to 

IDDZZ = 1 (default) standard treatment of FSR QCD corrections 

NPAR(22) = IAMW2 ZUFLAG( ^AMW2' , IAMW2) - incorporates the fermionic 
two-loop contributions to the prediction for the W boson mass [29,30,32]: 
IAMW2 = (default) no two-loop corrections to M^y 
IAMW2 = 1 with two-loop corrections to 

NPAR(23) = ISFSR ZUFLAG( ' SFSR' , ISFSR) - allows to switch the final 
state radiation 

ISFSR =-1 both, QED and QED®QCD final state radiation are excluded 
ISFSR = final state QED radiation is excluded, QED®QCDis included 
ISFSR = 1 (default) both, QED and QED®QCD final state radiation are 
included 



NPAR(24) = IDMWW ZUFLAG( 'DMWW , IDMWW) - Simulation of the theoreti- 
cal error on for AMT4 = 5,6: 
IDMWW = -1 minimal value for M^y 
IDMWW = (default) no shift on applied 
IDMWW = 1 maximal value for 

NPAR(25) = IDSWW ZUFLAG( 'DSWW' , IDSWW) - Simulation of the theoreti- 
cal error on sin^ 9^^^^ for AMT4 = 6: 
IDSWW = -1 minimal value for sin^ ^^g"* 
IDSWW = (default) no shift on sin^ 9^^^^ applied 
IDSWW = 1 maximal value for sin^ ^i?P' 
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A. 4 Calculation ofa{s). Function XF0TF3 

the running QED coupling at scale s is calculated with function XF0TF3 as 
follows: 

a 

a{s) = g . 

1 - — DREAL(XF0TF3(IALEM, IALE2, IHVP, IQCD, 1,DAL5H, -S)) 

4:71 

(A.2) 



B ZFITTER user guide 

This Appendix describes technical details of the ZFITTER package. Not all of 
them have been influenced by the updates over the years. For the convenience 
of the user, we nevertheless decided to give a complete overview and will repeat 
a substantial part of the material, which was already presented in Section 4.2 
of [2]. 

ZFITTER is coded in FORTRAN 77. Double-precision variables have been used 
throughout the program. The package consists of the following FORTRAN files: 

zf 630_aux.f 

zfbib6_40.f 

zfmai6_42.f 

zfusr6_42.f 

acol6_lp.f 

bcqcdl5_14.f 

bhang4_640.f 

bkqcdl5_14.f 

dizet6_42.f 

expif i6_30.f 

f unang6_30.f 

m2tcor5_ll.f 

pairho6_40.f 

APV_lib.f 
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The following routines are normally called in the initialization phase of pro- 
grams using the ZFITTER package in the order listed below: ZUINIT, ZUFLAG, 
ZUWEAK, ZUCUTS, ZUINFO. An example of different use is described in Section 
2.5 of [2]. 

B.l Subroutine ZmmT 

Subroutine ZUINIT is used to initialize variables with their default values. This 
routine must be called before any other ZFITTER routine. 



CALL ZUINIT 



B.2 Subroutine ZUFLAG 

Subroutine ZUFLAG is used to modify the default values of flags which control 
various ZFITTER options. 



CALL ZUFLAG (CHFLAG,IVALUE) 



Input Arguments: 

CHFLAG is the character identifier of a ZFITTER flag. 

IVALUE is the value of the flag. See Tab. C.2 for a list of the defaults. 

Possible combinations of CHFLAG and IVALUE are listed below: ^ 

In Tab. C.2 an overview over all flags used in DIZET and ZFITTER is given. 
The DIZET flags (vector NPAR(1 :25) in DIZET corresponds to NPARD(1 :25) in 
ZFITTER) are described in Subsection A. 

AFBC - Controls the calculation of the forward backward asymmetry for in- 
terfaces ZUTHSM, ZUXSA, ZUXSA2, and ZUXAFB: 

IVALUE = asymmetry calculation is inhibited (can speed up the program 
if asymmetries are not desired) 

^ It is worth noting that not for all flags the default value is necessarily the preferred 
value. A typical example is flag FINE, distinguishing two different treatments of FSR, 
which are relevant in different experimental setups. 



36 



IVALUE = 1 (default) both cross-section and asymmetry calculations are 
done 

AFMT - see NPAR(15) in subsection A.3 

ALEM - Controls the treatment of the running QED coupling a{s): 

IVALUE = or 2 Aa^^^ (Mj) must be supplied by the user as input to the 

DIZET package; using this input DIZET calculates ALQED = q;(M|) 
IVALUE = 1 or 3 Aa^^\M'^) and a(Mj) are calculated by the program 
using a parameterization IHVP. 

The scale of a (scale) is governed in addition by the flag CONV, see description 
below and the flowchart, Figure 6 of [2]. Values ALEM = 0,1 are accessible 
only at CONV = 0. Then for ALEM = 0,1 Q!(M|) and for ALEM = 2,3 a{s) are 
calculated. Values ALEM = 2,3 are accessible for CONV = 0,1,2. Then for CONV 
= a{s) and for CONV = 1,2 a{s') are calculated. Recommended values: ALEM 
= 2,3. 

ALE2 - see NPAR(19) in subsection A.3 
AMT4 - see NPAR(2) in subsection A.3 

ASCR - is a hidden flag that handles the treatment of contributions to ApB'- 
IVALUE = treatment as in versions 5 up to version 6.23 
IVALUE = 1 (default) new (and very old) treatment 

BARB - see NPAR(IO) in subsection A.3 

BORN - Controls calculation of QED and Born observables: 
IVALUE = (default) QED convoluted observables 

IVALUE = 1 electroweak observables corrected by Improved Born Approx- 
imation 

BOXD - Determines calculation of ZZ and WW box contributions, (see Section 
3.3.2 of [2]): 

IVALUE = no box contributions are calculated 

IVALUE = 1 (default) the boxes are calculated as additive separate contri- 
bution to the cross-section 
IVALUE = 2 box contributions are added to all four form factors 

CONV - Controls the energy scale of running a and EWRC, see Figure 6 of [2]: 
IVALUE = a{s) 

IVALUE = 1 (default) a{s') convoluted 

IVALUE = 2 both electroweak radiative correction and are convoluted 
CZAK - Treatment of CKHSS non-factorized corrections, [78], [79], see Figure 
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6 of [2]: 

IVALUE = without CKHSS corrections 
IVALUE = 1 (default) with CKHSS corrections everywhere 
IVALUE = 2 CKHSS corrections are taken into account only in POs, this 
option is used for tests only 

DIAG - Selects type of diagrams taken into account: 

IVALUE = —1 only Z- exchange diagrams are taken into account 
IVALUE = Z and 7 - exchange diagrams are taken into account 
IVALUE = 1 (default) Z and 7 exchange and Z7 interference are included 

EXPF - see NPAR(13) in subsection A.3 

EXPR - see NPAR(12) in subsection A.3 

ENUE - Treatment of the improved Born approximation for the process e^e~ — > 
Vf.Vf.. This option is only available via subroutine COSCUT in zfbib6_34.f 
and using zf EENN_34 . f . This flag is 'presently not supported, see also Section 
3. 

IVALUE = —1 s-channel only 
IVALUE = s- and t-channels 

IVALUE = 1 (default) s+t complete with s-t interference 
DMWW - see NPAR(24) in subsection A.3 
DSWW - see NPAR(25) in subsection A.3 

FBHO - treatment of second order corrections to angular distributions and ApB 
IVALUE = (default) treatment as in version 6.21 and before 
IVALUE = 1 modified treatment; photonic radiative corrections are improved 
and 0{a^) contributions from pairs are included in leading-log approxi- 
mation, see Section 2.1 

FINR - Controls the calculation of final-state radiation, 

IVALUE = —1 final-state QED and QCD correction are not applied; 
IVALUE = by s' cut, final-state QED correction is described with the 

factor 1 + 3a{s)/{47r)Qj 
IVALUE = 1 (default) M|j cut, includes complete treatment of final-state 

radiation with common soft-photon exponentiation 

F0T2 - Controls second-order leading log and next-to-leading log QED cor- 
rections: 

IVALUE — —1 no initial state radiation QED convolution at all 
IVALUE = complete a additive radiator 
IVALUE = 1 with logarithmic hard corrections 
IVALUE = 2 complete additive radiator 
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IVALUE = 3 (default) complete additive radiator 

IVALUE = 4 optional additive radiator for estimation of theoretical er- 
rors [89] 

IVALUE — 5 "pragmatic" LLA third order corrections in a factorized form [8] 

FSPP - correction due to final state radiation into pairs, see Section 2.4: 
IVALUE = (default) no final state state pair corrections 
IVALUE = 1 final state state pair contributions are implemented as additive 
corrections 

IVALUE = 2 final state state pair contributions are implemented as multi- 
plicative corrections 

FSRS - Final state radiation scale: 

IVALUE = a(0), preferred for tight cuts 

IVALUE = 1 (default) a(s), preferred for loose cuts 

FT JR - Treatment of FT JR corrections [77] : 
IVALUE = without FTJR corrections 
IVALUE = 1 (default) with FTJR corrections everywhere 
IVALUE = 2 FTJR corrections are taken into account only in POs, the op- 
tion is used for tests only 

FUNA - implementation of higher order photonic corrections to the angular 
distribution, see Section 2.2: 

IVALUE — (default) no higher order corrections 

IVALUE = 1 higher order photonic LLA corrections are included 

GAMS - Controls the s dependence of Qz^ the Z-width function, see Section 
3.2 of [2]: 

IVALUE = forces Qz to be constant 

IVALUE = 1 (default) allows Qz to vary as a function of s [90]. 
GFER - see NPAR(20) in subsection A.3 
HIGS - see NPAR(14) in subsection A.3 
HIG2 - see NPAR(18) in subsection A.3 

INTF - Determines if the 0{q) initial-final state QED interference (IFI) is 
calculated; see Section 2.3: 
IVALUE = the interference term is ignored 
IVALUE = 1 (default) with IFl in the O {a) 

IVALUE = 2 with one-loop IFI corrections and corresponding higher order 
effects from the exponentiation 

IPFC - Pair flavour content for the pair production corrections: 
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IVALUE = 1 only electron pairs 

I VALUE = 2 only muon pairs 

IVALUE = 3 only tau-lepton pairs 

IVALUE = 4 only hadron pairs 

IVALUE = 5 (default) all channels summed 

IVALUE = 6 leptonic pairs (without hadrons) 

IPSC - Pair production singlet-channel contributions (works with ISPP = 2): 
IVALUE = (default) only non-singlet pairs 
IVALUE = 1 LLA singlet pairs according to [91] 
IVALUE = 2 complete O(a^) singlet pairs, ibid 
IVALUE = 3 singlet pairs up to order (aL)^, ibid 

IPTO - Third (and higher) order pair production contributions [12] (works 
with ISPP = 2): 

IVALUE = —1 allows to calculate the pure virtual pair contributions sepa- 
rately 

IVALUE = only O(a^) contributions 
IVALUE = 1 0{a^) pairs 

IVALUE = 2 some "non- standard" 0{a^) LLA pairs added 
IVALUE = 3 (default) O(a^) LLA electron pairs added 

ISPP - Treatment of ISR pairs: 

IVALUE = —1 pairs are treated with a "fudge" factor as in versions up to 
v.5.14 

IVALUE = without ISR pairs 

IVALUE = 1 with ISR pairs, [92] with a re-weighting 
IVALUE = 2 (default) with ISR pairs according to [12] 
IVALUE = 3 with ISR pairs according to [93] 

IVALUE = 4 with ISR pairs [93] with extended treatment of hadron pair 
production 

MISC - Controls the treatment of scaling of p in the Model Independent ap- 
proach, see discussion in Subsection 6.1: 
IVALUE = (default) non-scaled p's are used, AROTFZ-array 
IVALUE = 1 scaled p's, absorbing imaginary parts, are used, ARROFZ-array 

MISD - Controls the s dependence of the Model Independent approach 
IVALUE = fixed s = M| in EWRC, old treatment 

IVALUE = 1 (default) ensures equal numbers from all interfaces and for all 
partial channels but INDF=10 for a large range of ^/s and for INDF=10 up 
to 100 GeV 

PART - Controls the calculation of various parts of Bhabha scattering: 

IVALUE = (default) calculation of full Bhabha cross-section and asym- 
metry 
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I VALUE = 1 only s channel 
I VALUE = 2 only t channel 
I VALUE = 3 only s — t interference 

POWR - Controls inclusion of final-state fermion masses in kinematical factors, 
see Figure 6 of [2]. It acts differently for quarks and leptons. For leptons: 
IVALUE = final state lepton masses are set equal to zero 
IVALUE = 1 (default) final state lepton masses are retained in all kinemat- 
ical factors 

For quarks it is active only for FINE = -1 and then: 
IVALUE = final state quark masses are set equal to zero 
IVALUE = 1 (default) final state quark masses are set to their running val- 
ues (that is, for cc and bb channels) and retained in all kinematical factors 

PREC - is an integer number which any precision governing any numerical 
integration is divided by, increasing thereby the numerical precision of com- 
putation: 

IVALUE = 10 (default) 

IVALUE = 1 — 99 in some cases when some numerical instability while run- 
ning v.5.10 was registered, it was sufficient to use PREC = 3, in some other 
cases (e.g. with Pr) only PREC = 30 solved the instability 

PRNT - Controls ZUWEAK printing: 

IVALUE = (default) printing by subroutine ZUWEAK is suppressed 
IVALUE = 1 each call to ZUWEAK produces some output 

QCDC - see NPAR(3) in subsection A.3 

SCAL - see NPAR(9) in subsection A.3 

SCRE - see NPAR(6) in subsection A.3 

SFSR - see NPAR(23) in subsection A.3 

TUPV - simulates theoretical uncertainties in APV, see Sect. 4: 
IVALUE = 1 (default) 
IVALUE = 2,3 for variation, see Sect. 4 

VPOL - see NPAR(l) in subsection A.3 

WEAK - Determines if the weak loop calculations are to be performed 

IVALUE = —1 is only valid for INDF=— 1 and represents the the options 
WEAK=1,2 without electroweak corrections for W boson exchange in the 
t-channel 

IVALUE = no weak loop corrections to the cross-sections are calculated 
and weak parameters are forced to their Born values, i.e. pe/ = i^ej,ef = 1 
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IVALUE = 1 (default) weak loop corrections to the cross-sections are cal- 
culated 

IVALUE = 2 weak loop corrections are calculated but some higher order 
corrections that do not propagate via DIZET are swiched off (ADDIME, 
ADDIMF; these corrections are small at LEP2 and they are not included 
by using DIZET with other codes 

B.3 Subroutine ZUWEAK 

Subroutine ZUWEAK is used to perform the weak sector calculations. These are 
done internally with DIZET, see Section A. The routine calculates a number 
of important electroweak parameters which are stored in common blocks for 
later use (see Section 6.1). If any ZFITTER flag has to be modified this must 
be done before calling ZUWEAK. 



CALL ZUWEAK (ZMASS , TMASS , HMASS , DAL5H , ALFAS) 



Input Arguments: 

ZMASS is the Z mass in GeV. 
TMASS is the top quark mass mj in GeV, [10-400]. 
HMASS is the Higgs mass in GeV, [10-1000]. 
DAL5H is the value of Aal^\Ml). 

ALFAS is the value of the strong coupling constant at = M| (see factors 
QCDCOR in Tab. A.l). 

Computing time may be saved by performing weak sector calculations only 
once during the initialization of the ZFITTER package. This is possible since 
weak parameters are nearly independent of s near the Z peak, e.g. ~ In s/M|. 
However, the incredible precision of LEPl data forced us to give up this option, 
see description of flag MISD. 

B.4 Subroutine ZUCUTS 

Subroutine ZUCUTS is used to define kinematic and geometric cuts for each 
fermion channel: it selects the appropriate QED calculational chain. 



CALL ZUCUTS ( INDF , ICUT , ACOL , EMIN , S_PR , ANGO , ANGl , SIPP) 
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Input Arguments: 



INDF is the fermion index (see Tab. C.l and Figure 5 of [2]). 
I CUT controls the kinds of cuts (chain) to be used. 

ICUT = —1: (defauh) allows for an s' cut (a cut on Mj^, the fermion and 

antifermion invariant mass); the fastest branch based on [94] 
ICUT = 0: not recommended; branch is known to contain bugs. It allows 
for a cut on the acollinearity ACOL of the // pair, on the minimum energy 
EMIN of both fermion and antifermion, and for a geometrical acceptance 
cut [95] 6 

ICUT = 1: s' or cuts and geometrical acceptance cut, based on [99] 
ICUT = 2: new branch, replaces ICUT = for realistic cuts ACOL and EMIN, 
based on [96] 

ICUT = 3: the same branch, using ACOL cut and EMIN cut but also with 
possibility to impose an additional acceptance cut [98] 
ACOL is the maximum acollinearity angle ^"^'"^^ of the // pair in degrees (ICUT 
= 0,2,3). 

EMIN is the minimum energy E^"^ of the fermion and antifermion in GeV 
(ICUT = 0,2,3). 

S_PR is the minimum allowed invariant // mass M^j in GeV (ICUT = -1,1) 
or, with some approximations, the minimum allowed invariant mass of the 
propagator after ISR^ 

ANGO (default = 0°) is the minimum polar angle in degrees of the final-state 
antifermion. 

ANGl (default = 180°) is the maximum polar angle in degrees of the final- 
state antifermion. 

SIPP is a parameter for cuts on the invariant mass of secondary pairs for 
FSPP=1,2 (see subsection 2.4). In older versions before v.6.30 the parameter 
SIPP governed the calulation of corrections due to initial state pairs and 
was recommended to be chosen equal to s'. 

B.5 Subroutine ZUINFO 

Subroutine ZUINFO prints the values of ZFITTER flags and cuts. 



CALL ZUINFO (MODE) 



^ As was shown recently [96-98], the old results of [95] contained bugs which oc- 
casionally didn't show up in comparisons as e.g. in [39]. The option is retained for 
back-compatibility with older versions only. 

^ The invariant mass of the propagator is not an observable quantity unless specific 
assumptions on ISR and FSR are made. 
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Input Argument: 



MODE controls the printing of ZFITTER flag and cut values. 
MODE = 0: Prints all flag values. 
MODE = 1: Prints all cut values. 



C Interface routines of ZFITTER 

This Appendix describes technical details of the interface routines of ZFITTER. 
Not all of them have been influenced by the updates over the years. For the 
convenience of the user, we nevertheless decided to give a complete overview 
and will repeat a substantial part of the material, which was already presented 
in Section 4.3 of [2]. 

Note that subroutine ZUWEAK must be called prior to the interfaces. As a 
consequence, the flags used in ZUWEAK may influence the calculation of cross- 
sections and asymmetries in the interfaces described now. 

All subroutines need the following input arguments: 

INDF is the fermion index (see Tab. C.l). 
SQRS is the centre-of-mass energy ^/s in GeV. 
ZMASS is the Z mass Mz in GeV. 

We just mention that the interface using an S-matrix inspired language is 
reahzed with the SMATASY package [100-103]. 

C.l SubrouUne ZUTESn 

Subroutine ZUTHSM is used to calculate Standard Model cross-sections and 
forward-backward asymmetries. 



CALL ZUTHSM ( INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , ALFAS , XS* , AFB* ) 



Input Arguments: 

TMASS is the top quark mass rrit in GeV, [10-400]. 
HMASS is the Higgs mass in GeV, [10-1000]. 
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DAL5H is the value of Aajf^ (M|). 

ALFAS is the value of the strong coupling constant at = M| (see also 
flag QCDC and factors QCDCOR). 

Output Arguments ^ : 

XS is the total cross-section ctt in nb. 

AFB is the forward-backward asymmetry A^^. 

Output Internal Flag: 
INTRF=1 



C.2 Subroutine ZUATSM 

Subroutine ZUATSM is used to calculate differential cross-sections, da/dcos6, 
in the Standard Model. 



CALL ZUATSM ( INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , ALFAS , CSA* , DXS* ) 



Input Arguments: 

TMASS is the top quark mass rrit in GeV, [10-400]. 
HMASS is the Higgs mass in GeV, [10-1000]. 

ALQED is the value of the running electromagnetic coupling constant. 
ALFAS is the value of the strong coupling constant ag at = M| (see factors 
QCDCOR). 

CSA is the cosine of the scattering angle. 
Output Arguments: 

DXS is the theoretical differential cross- sect ion. 
Output Internal Flag: 

INTRF=1 



^ An asterisk (*) following an argument in a calling sequence is used to denote an 
output argument. 
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C.3 Subroutine ZUTPSM 



Subroutine ZUTPSM is used to calculate the tau polarization and tau polariza- 
tion asymmetry in the Standard Model. 



CALL ZUTPSM (SQRS , ZMASS , TMASS , HMASS , DAL5H , ALFAS , TAUPOL* , TAUAFB*) 



Input Arguments: 

HMASS is the Higgs mass in GeV, [10-1000]. 
DAL5H is the value of Aa'j^\Ml). 

ALFAS is the value of the strong coupling constant ag at = M| (see factors 
QCDCOR). 

Output Arguments: 

TAUPOL is the tau polarization Apoi of Equation 3.312 of [2]. 
TAUAFB is the tau polarization forward-backward asymmetry A^°g as defined 
in Equation 3.313 of [2]. 

Output Internal Flag: 
INTRF=1 

C.4 Subroutine ZULRSM 

Subroutine ZULRSM is used to calculate the left-right asymmetry in the Stan- 
dard Model. 



CALL ZULRSM ( INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , ALFAS , POL , XSPL* , XSMI * ) 



Input Arguments: 

TMASS is the top quark mass rut in GeV, [40-300]. 
HMASS is the Higgs mass in GeV, [10-1000]. 
DAL5H is the value of Aa'i\Ml). 
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ALFAS is the value of the strong couphng constant as at 
flag ALST). 

POL is the degree of longitudinal polarization of electrons. 
Output Arguments: 

XSPL is the cross-section for POL > 
XSMI is the cross-section for POL < 

Output Internal Flag: 
INTRF=1 



C.5 Subroutine ZUXSA 

Subroutine ZUXSA is used to calculate cross-section and forward-backward 
asymmetry as described in Section 3.7.5 of [2] as functions of y^, M^, T z- 



CALL ZUXSA ( INDF , SQRS , ZMASS , GAMZO , MODE , GVE , XE , GVF , XF , XS* , AFB* ) 



Input Arguments: 

GAMZO is the total Z width in GeV. 
MODE determines which weak couplings are used: 

MODE = 0: XE (XF) are effective axial-vector couplings Cej for electrons 

(final-state fermions). 
MODE = 1: XE (XF) are the effective weak neutral-current amplitude nor- 
malizations pej for electrons (final-state fermions). 
GVE is the effective vector coupling for electrons Qe- 

XE is the effective axial-vector coupling a^. or weak neutral-current amplitude 
normalization for electrons (see MODE). 

GVF is the effective vector coupling for the final-state fermions pf. 

XF is the effective axial- vector coupling a/ or the weak neutral-current ampli- 
tude normalization pf for the final-state fermions (see MODE). 

Output Arguments: 

XS is the cross-section in nb. 

AFB is the forward-backward asymmetry A^^. 



= M| (see also 
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Output Internal Flag: 



INTRF=3 



C.6 Subroutine ZUXSA2 

Subroutine ZUXSA2 is used to calculate lepton cross-section and forward-backward 
asymmetry as functions of -y/i, M^j, F^, and of the weak couplings assuming 
lepton universality. This routine is similar to ZUXSA except that the couphngs 
are squared. 



CALL ZUXSA2 ( INDF , SQRS , ZMASS , GAMZO , MODE , GV2 , X2 , XS* , AFB*) 



Input Arguments: 

GAMZO is the total Z width in GeV. 

MODE determines which weak couplings are used: 

MODE = 0: X2 is the square of the effective axial-vector coupling ai for lep- 
tons. 

MODE = 1: X2 is the square of the effective neutral-current amplitude nor- 
malization pi for leptons. 
GV2 is the square of the effective vector coupling gi for leptons. 
X2 is the square of the effective axial-vector coupling ai or neutral-current 

amplitude normalization pi for leptons (see MODE). 

Output Arguments: 

XS is the cross-section (Tt in nb. 

AFB is the forward-backward asymmetry A^g. 

Output Internal Flag: 

INTRF=4 



C. 7 Subroutine ZUTAU 

Subroutine ZUTAU is used to calculate the polarization as a function of ^i, 
M^, F^, and the weak couplings. 
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CALL ZUTAU (SQRS , ZMASS , GAMZO , MODE , GVE , XE , GVF , XF , TAUPOL* , TAUAFB*) 



Input Arguments: 

GAMZO is the total Z width in GeV. 

MODE determines which weak couphngs are used: 

MODE = 0: XE (XF) is the effective axial-vector coupling ciej for electrons 

(final-state fermions). 
MODE = 1: XE (XF) is the effective weak neutral-current amplitude normal- 
ization pe,/ for electrons (final-state fermions). 
GVE is the effective vector coupling for electrons ge ■ 

XE is the effective axial-vector coupling ae or weak neutral-current amplitude 

normalization pe for electrons (see MODE). 
GVF is the effective vector coupling for the final-state fermions . 
XF is the effective axial- vector coupling or weak neutral-current amplitude 

normalization p/ for the final-state fermions (see MODE). 

Output Arguments: 

TAUPOL is the r-lepton polarization A7- 

TAUAFB is the forward-backward asymmetry for polarized r-leptons 

Output Internal Flag: 

INTRF=3 



C.8 Subroutine Z\]XSEC 

Subroutine ZUXSEC is used to calculate the cross section as a function of y/s, 
Mz, Tz, Te and Tf. 



CALL ZUXSEC ( I NDF , S QRS , ZMASS , GAMZO , GAMEE , GAMFF , XS * ) 



Input Arguments: 

GAMZO is the total Z width Tz in GeV. 
GAMEE is the partial Z decay width Fg in GeV. 

GAMFF is the partial Z decay width Tf in GeV; if INDF=10, GAMFF=F,,. 
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Output Internal Flag: 



INTRF=2 

C.9 Subroutine ZUXAFB 

Subroutine ZUXAFB is used to calculate the cross section as a function of a/s, 
Mz, Tz, Te and Tf. 

CALL ZUXAFB ( INDF , SQRS , ZMASS , GAMZO , PFOUR , PVAE2 , PVAF2 , XS* , AFB* ) 
Input Arguments: 

GAMZO is the total Z width in GeV. 

PFOUR is the product of vector and axial- vector couplings geaeQfaf. 
PVAE2 is gj + al- 
PVAF2 is gj + aj. 

Output Argument: 

XS is the cross-section cr^ in nb. 

AFB is the forward-backward asymmetry A^^. 

Output Internal Flag: 
INTRF=5 

C.IO Subroutine ZUALR 

Subroutine ZUALR is reserved for the fit of ^4^^. 

CALL ZUALR (SQRS , ZMASS , GAMZO , MODE , GVE , XE , GVF , XF , TAUPOL* , TAUAFB*) 

Output Internal Flag: 
INTRF=6 
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C.ll Subroutine ZVVEkK 



The strength of the Wtb vertex, \ Vtb\, can be changed by using the subroutine 
ZVWEAK instead of ZUWEAK; see Section 1. 



CALL ZVWEAK ( ZMASS , TMASS , HMASS , DAL5H , V_TB , ALFAS ) 



Input Arguments: 

V_TB is the value of \Vtb\', the Standard Model calculation of ZUWEAK corre- 
sponds to \Vtb\ = 1. 

ZVWEAK has to be complemented with new interface routines, all now with the 
the additional input argument V_TB: 

ZVTHSM ( INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , V_TB , ALFAS , XS* , AFB*) 
ZVTPSM (SQRS , ZMASS , TMASS , HMASS , DAL5H , V_TB , ALFAS , TAUPOL* , TAUAFB* ) 
ZVLRSM (INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , V_TB , ALFAS , POL , XSPL* , XSMI*) 
ZVATSM (INDF , SQRS , ZMASS , TMASS , HMASS , DAL5H , V_TB , ALFAS , CSA , DXS*) 

These subroutines replace ZUTHSM, ZUTPSM, ZULRSM and ZUATSM. 

C.12 Subroutine ZU_kPY 

With subroutine ZU_APV measurements of the weak charge, Qw^ (see Section 4) 
can be included into the global tests of the Standard Model. 



CALL ZU_APV (ZMASS , TMASS , HMASS , SIN2TW , UMASS , DMASS , ClU* , CID* , C2U* , C2D*) 



Input Arguments: 

SIN2TW is the sin of the weak mixing angle. 
UMASS is the u-quark mass (constituent). 
DMASS is the d-quark mass (constituent). 

Output Argument: 

ClU - coupling parameter, Ciu = 2aeW„ 
CID - coupling parameter. Cm = 2aeVd 
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C2U - coupling parameter, C2« = 2fea„ 
C2D - coupling parameter, C2d = 2We'3d 
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Final- 
state 
fermions 


vv e^e 




r+r- 


uu dd cc 


ss tt bb 


hadrons 


Bhabha 


INDF 


1 


2 


3 


4 5 6 


7 8 9 


10 


11 



Table C.l 

Indices for the selection of final states. Note that INDF=0 returns values for one 
neutrino species, while INDF=10 returns values for the five-flavour inclusive (udsch) 
hadronic channel. Also note that INDF = 0,1 includes only s-channel calculations 
while INDF = 8 always returns zero. 
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I 


' FLAG ' 


name m 
programs 


Position in DIZET 
NPARd :25) 


Position in ZFITTER 
NPARd: 30) 


default 
value 


1 


AFBC 


lAFB 




13 


1 


2 


SCAL 


ISCAL 


9 


15 





3 


SCRE 


ISCRE 


6 







4 


AMT4 


IAMT4 


2 


16 


4 


5 


BORN 


IBQRN 




14 





6 


BOXD 


IBOX 




4 


1 


7 


CONV 








1 


8 


FINR 


IFINAL 




9 


1 


9 


F0T2 


IPH0T2 




10 


3 


10 


GAMS 






5 


1 


11 


DIAG 






7 


1 


12 


INTF 


INTERF 




8 


1 


13 


BARB 


IBARB 


10 




2 


14 


PART 


IPART 




17 





15 


POWR 








1 


16 


PRNT 











17 


ALEM 


lALEM 


7 


20 


3 


18 


QCDC 


IQCD 


3 


3 


3 


19 


VPOL 


IHVP 


1 


2 


1 


20 


WEAK 


IWEAK 




1 


1 


21 


FTJR 


IFTJR 


11 




1 


22 


EXPR 


IFACR 


12 







23 


EXPF 


IFACT 


13 


19 





24 


HIGS 


IHIGS 


14 







25 


AFMT 


lAFMT 


15 




3 


26 


CZAK 


ICZAK 


17 




1 


27 


PREC 


NPREC 






10 


28 


HIG2 


IHIG2 


18 







29 


ALE2 


IALE2 


19 


21 


3 


30 


GFER 


IGFER 


20 




2 


31 


ISPP 


ISRPPR 






2 


32 


FSRS 








1 


33 


MISC 


IMISC 









34 


MISD 


IMISD 






1 


35 


IPFC 








5 


36 


IPSC 











37 


IPTO 








3 


38 


FBHO 











39 


FSPP 


IFSPPR 









40 


FUNA 


IFUNAN 









41 


ASCR 








1 


42 


SFSR 


ISFSR 


23 




1 


43 


ENUE 








1 


44 


TUPV 


ITUPV 






1 


45 


DMWW 


IDMWW 


24 







46 


DSWW 


IDSWW 


25 











IMOMS 


4 




1 






IMASS 


5 











IMASK 


8 











lEWLC 


16 




1 






IDDZZ 


21 




1 






IAMW2 


22 








Table C.2 

Flag settings for ZFITTER and DIZET; the flags are listed in the order of vector 
IFLAGS. The corresponding names used internally in the programs, the position of 
the flags in vector NPAR(1 : 25) o/DIZET fca/Zed NPARD mZFITTERj and NPAR(1 : 30) 
of ZFITTER and the default values are given 



54 



References 



[1] D. Bardin, P. Christova, M. Jack, L. Kalinovskaya, A. Olshevski, S. Riemann, 
and T. Riemann, Fortran program ZFITTER version 6.21 (26 July 1999), 
obtainable from webpage 

[http: / /www-zeuthen. desy.de/tlieory/research/zfitter/ 



[2] D. Bardin, P. Christova, M. Jack, L. Kalinovskaya, A. Olshevski, S. Riemann, 
and T. Riemann, Comput. Phys. Commun. 133 (2001) 229-395, 
hep-ph/9908433| 



[3] ZFITTER support group, A. B. Arbuzov, M. Awramik, M. Czakon, A. 
Freitas, M. W. Griinewald, K. Monig, S. Riemann, T. Riemann, webpage 
[http : / / www-zeuthen . desy. de / theory / r esearch/zfitter / 1 



[4] D. Bardin, M. Griinewald, and G. Passarino, |hep-ph/9902452 



[5] D. Bardin and G. Passarino, The standard model in the making: Precision 
study of the electroweak interactions. Clarendon, Oxford, UK, 1999. 



[6] Two Fermion Working Group Collaboration, M. Kobel et al., hep-ph/0007180 



[7] LEP Electroweak Working Group, webpage http: / /lepewwg.web.cern.ch/ 
LEPEWWG/. 

[8] M. Skrzypek, Acta Phys. Polon. B23 (1992) 135. 

[9] A. Arbuzov, Phys. Lett. B470 (1999) 252-258, plep-ph/990836l1 

[10] D. Bardin, M. Bilenky, A. Chizhov, A. Sazonov, Y. Sedykh, T. Riemann, and 
M. Sachwitz, Phys. Lett. B229 (1989) 405. 

[11] W. Beenakker, F. Berends, and W. van Neerven, "Applications of 
renormalization group methods to radiative corrections", in Proc. of the 
Int. Workshop on Radiative Corrections for e'^e~ Collisions, Schlofi Ringberg, 
Tegernsee, Germany, April 1989 (J. Kiihn, ed.), p. 3, Springer, Berlin, 1989. 



[12] A. Arbuzov, JHEP 07 (2001) 043 and hep-ph/9907500 



[13] M. Greco and A. Grillo, Nuovo Cim. Lett. 15 (1976) 174. 

[14] S. Jada ch, B. Ward, and Z. Was, Phys. Rev. D63 (2001) 113009, 
jhep-ph/0 006359, 

[15] A. Hoang, J. Kiihn, and T. Teubner, Nucl. Phys. B455 (1995) 3-20, 
hep-ph/9507255, 

[16] S. Jadach, B. Ward, and Z. Was, Comput. Phys. Commun. 130 (2000) 260- 
325, hep-ph /9912214, 

[17] D. Bardin, S. Jadach, T. Riemann, and Z. Was, Eur. Phys. J. C24 (2002) 
373-383, hep-ph/0110371 . 



55 



[18] D. Bardin, C. Burdik, P. Khristova, and T. Riemann, Z. Phys. C44 (1989) 
149. 



[19] D. Bardin, C. Burdik, P. Khristova, and T. Riemann, "Study of electroweak 
radiative corrections to deep inelastic scattering at HERA", subm. to: 
Workshop on Physics at HERA, DESY, Hamburg, Oct 12-14, 1987, Dubna 
preprint JINR E2-87-595 (1987). 

[20] D. Bardin, P. Christova, L. Kahnovskaya, and G. Passarino, Eur. Phys. J. 
C22 (2001) 99-104, hep-ph/0102233. 



[21] A. Andonov et ah, Phys. Part. Nucl. 34 (2003) 577-618, pIep-ph/0207156 

[22] W. Marciano and A. Sirhn, Phys. Rev. D27 (1983) 552. 

[23] G. Degrassi, S. Fanchiotti, F. Ferugho, B. Gambino, and A. Vicini, "Two- 
loop electroweak top corrections: are they under control?" , in Reports of the 
Working Group on Precision Calculations for the Z Resonance, report CERN 
95-03 (1995) (D . Bardin, W. Hollik, and G. Passarino, eds.), pp. 163-174, 
hep-ph/9412380| 



[24] G. Degrassi, S. Fanchiotti, F. Feruglio, P. Gambino, and A. Vicini, "Two loop 
corrections to the heavy top limit in the p parameter" , New York Univ. preprint 
NYU-TH-01-04-95 (1995). 

[25] G. Degrassi, F. Feruglio, A. Vicini, S. Fanchiotti, and P. Gambino, "Two 
loop corrections for electroweak processes", in Proc. of the 30*^* Rencontre 
de Moriond, Les Arcs, France, March 1995 (J. Tran Than Van, ed.), vol. 1, 
pp. 77-86, Editions Frontieres, 1995, hep-ph/9507286,. 

[26] G. Degrassi, P. Gambino, and A. Vicini, Phys. Lett. B383 (1996) 219-226, 
hep-ph/9603374 , 

[27] G. Degrassi, Fortran program m2tcor (Oct 1996). 

[28] G. Degrassi and P. Gambino, Nucl. Phys. B567 (2000) 3-31, |hep-ph/9905472| 
[29] A. Freitas, W. Holhk, W. Walter, and G. Weiglein, Phys. Lett. B495 (2000) 



338-346, hep-ph/0007091, E: B570 (2003) 260 



[30] A. Freitas, S. Heinemeyer, W. Hollik, W. Walter, and G. Weiglein, Nucl. Phys. 
Proc. Suppl. 89 (2000) 82-87, |hep-ph/0007"l29l 

[31] A. Freitas, S. Heinemeyer, W. Hollik, W. Walter, and G. Weiglein, 
|hep-ph/0101260| 

[32] A. Freitas, W. Hollik, W. Walter, and G. Weiglein, Nucl. Phys. B632 (2002) 



189-218, E: B666 (2003) 305-307, |hep-ph/0202131 



[33] M. Awramik and M. Czakon, Phys. Lett. B568 (2003) 48-54, plep^h /0305248| 

[34] M. Awramik and M. Czakon, Phys. Rev. Lett. 89 (2002) 241801, 
.hep-ph/0208113, 



56 



[35] M. Awramik and M. Czakon, Nucl. Phys. Proc. Suppl. 116 (2003) 238-242, 
|hep-ph/0211041| 

[36] M. Awramik, M. Czakon, A. Onishchenko, and O. Veretin, Phys. Rev. D68 



(2003) 053004, hep-ph/0209084 



[37] M. Awramik, M. Czakon, A. Freitas, and G. Weiglein, Phys. Rev. D69 (2004) 
053006, [hep-ph/0311148t 

[38] A. Onishchenko and O. Veretin, Phys. Lett. B551 (2003) 111-114, 
|hep-ph/ 0209010t 

[39] D. Bardin, G. Passarino, and W. Holhk (eds.), "Reports of the working group 
on precision calculations for the Z resonance", report CERN 95-03 (1995). 

[40] A. Djouadi and C. Verzegnassi, Phys. Lett. B195 (1987) 265. 

[41] A. Djouadi, Nuovo Cim. AlOO (1988) 357. 

[42] B. A. Kniehl, Nucl. Phys. B347 (1990) 86-104. 

[43] F. Halzen and B. A. Kniehl, Nucl. Phys. B353 (1991) 567-590. 

[44] B. A. Kniehl and A. Sirlin, Nucl. Phys. B371 (1992) 141-148. 

[45] B. A. Kniehl and A. Sirhn, Phys. Rev. D47 (1993) 883-893. 

[46] A. Djouadi and P. Gambino, Phys. Rev. D49 (1994) 3499-3511, E: D53 (1996) 



4111, 



hep-ph/9309298 



[47] L. Avdeev, J. Fleischer, S. Mikhailov, an d O. Tar asov, Phys. Lett. B336 
(1994) 560-566, E: ibid., B349 (1995) 597, |hep-ph/9 406363, 

[48] K. Chet yrkin, J. Kiihn, and M. Steinhauser, Phys. Lett. B351 (1995) 331-338, 
l| hep-ph /950229l( 

[49] K. Chetyr kin, J. Kuhn, and M. Steinhauser, Phys. Rev. Lett. 75 (1995) 
3394-3397, |hep-ph/9504'4T3l 

[50] K. G. C hetyrkin, J. H. Ki ihn, and M. Steinhauser, Nucl. Phys. B482 (1996) 
|hep-ph/9606230| 



213-240, 



[51] M. Faisst, J. H. Kiihn, T. Seidensticker, and O. Veretin, Nucl. Phys. B665 
(2003) 649-662, |hep-ph/0302275| 

[52] M. Awramik, M. Czakon, A. Freitas, and G. Weiglein, Phys. Rev. Lett. 93 



(2004) 201805, hep-ph/0407317 



[53] M. Awramik, M. Czakon, A. Freitas, and G. Weiglein, Nucl. Phys. Proc. 
Suppl. 135 (2004) 119-123, plep-ph/0408207l 

[54] M. Awramik, M. Czakon, A. Freitas, and G. Weiglein, | hep-ph/0409T42l 

[55] R. Barbieri, M. Beccaria, P. Ciafaloni, G. Curci, and A. Vicere, Phys. Lett. 
B288 (1992) 95-98, E: B312 (1993) 511, pIep-ph/9205238 , 



57 



[56] R. Barbieri, M. Beccaria, P. Ciafaloni, G. Curci, and A. Vicere, Nucl. Phys. 
B409 (1993) 105-127. 

[57] J. Fleischer, O. Tarasov, and F. Jegerlehner, Phys. Lett. B319 (1993) 249-256. 

[58] J. Fleischer, O. Tarasov, and F. Jegerlehner, Phys. Rev. D51 (1995) 3820- 
3837. 

[59] G. Degrassi, P. Gambino, and A. Sirlin, Phys. Lett. B394 (1997) 188-194, 
|hep-ph/961 1363| 



[60] A. Freitas and K. Monig, Eur. Phys. J. C40 (2005) 493, |hep-ph/0411304 



[61] D. Bardin, P. Christova, and O. Fedorenko, "On the lowest order electroweak 
corrections to spin-^ fermion scattering (I). The one loop diagrammar" , Nucl. 
Phys. B175 (1980) 435. 

[62] D. Bardin, P. Christova, and O. Fedorenko, "On the lowest order electroweak 
corrections to spin-^ fermion scattering (II). The one-loop amplitudes", Nucl. 
Phys. B197 (1982) 1. 

[63] A. B. Arbuzov, D. Bardin, and A. Leike, Mod. Phys. Lett. A7 (1992) 2029- 
2038, E: ibid. A9 (1994) 1515. 

[64] D. Bardin, S. Riemann, and T. Riemann, Z. Phys. C32 (1986) 121. 

[65] S. Eidelman and F. Jegerlehner, Z. Phys. C67 (1995) 585-602, 
.hep-ph/9502298, 

[66] F. Jegerlehner, "Renormalizing the standard model" , in Testing the Standard 
Model: Proceedings of the Theoretical Advanced Study Institute in Elementary 
Particle Physics (TASI), Boulder, CO, Jun 3-29, 1990 (M. Cvetic and 
P. Langacker, eds.), pp. 476-590, World Scientific, Teaneck, N.J., 1991, 
[preprint PSl-PR-9r08l 

[67] F. Jegerlehner, Prog. Part. Nucl. Phys. 27 (1991) 1-76. 

[68] H. Burkhardt, F. Jegerlehner, G. Penso, and C. Verzegnassi, Z. Phys. C43 
(1989) 497. 

[69] M. Consoli, W. Hollik, and F. Jegerlehner, "Electroweak radiative corrections 
for Z physics", in Proc. of Workshop on Z Physics at LEP, Geneva, 
Switzerland, Feb 20-21 and May 8-9, 1989, report CERN 89-08 (1989) 
(G. Altarelli, R. Kleiss, and C. Verzegnassi, eds.), vol. 1, p. 7. 

[70] S. Fanchiotti and A. Sirlin, "Higher order contributions to Ar", in M.A.B. 
Beg Memorial Volume, Proc. Conf. on High Energy Physics and Cosmology, 
Islamabad, 1990 (A. Ah and P. Hoodbhoy, eds.), pp. 58-69, World Scientific, 
Singapore, 1991. 

[71] D. Bardin and A. Chizhov, "On the 0(a 

emO^s) corrections to electroweak 
observables" , in Proc. Int. Topical Seminar on Physics of e^e~ Interactions at 
LEP energies, 15-16 Nov 1988, JINR Dubna, USSR, JINR preprint E2-89-525 
(1989) (D. Bardin et al., eds.), pp. 42-48, Dubna prepr. E2-89-525, 



58 



[72] B. A. Kniehl, "Estimation of higher-order QCD effects on electroweak 
parameters", in Reports of the Working Group on Precision Calculations 
for the Z Resonance, report CE RN 95-03 (1995) (D. Bardin, W. Hollik, and 
G. Passarino, eds.), pp. 299-312, |he "p-ph/9503396l 

[73] A. Sirlin, "On the QCD corrections to in Reports of the Working Group 
on Precision Calculations for the Z Resonance, report CERN 95-03 (1995) 
(D. Bardin, W. Hollik, and G. Passarino, eds.), pp. 285-298, |hep^/9503396| 

[74] J. J. van der Bij and F. Hoogeveen, Nucl. Phys. B283 (1987) 477. 

[75] R. Barbieri, P. Ciafaloni, and A. Strumia, Phys. Lett. B317 (1993) 381. 

[76] R. Barbieri, private communication. 

[77] J. Fleischer, O. Tarasov, F. Jegerlehner, and P. Raczka, Phys. Lett. B293 
(1992) 437-444. 

[78] A. Czarnecki an d J. H. Kiihn, Phys. Rev. Lett. 77 (1996) 3955-3958, 
|hep-ph/9 608366f 

[79] R. Harlander, T. Seidensticker, and M. Steinhauser, Phys. Lett. B426 (1998) 
|hep-ph/97l'2228l 



125-132, 



[80] J. van der Bij and M. Veltman, Nucl. Phys. B231 (1984) 205. 

[81] J. J. van der Bij, Nucl. Phys. B248 (1984) 141. 

[82] G. Kallen and A.Sabry, K. Dan. Vidensk. Selsk. Mat.-Fys. Medd. 17 (1955) 
29. 

[83] M. Steinhauser, Phys. Lett. B429 (1998) 158-161, |hep-ph/9803313] 

[84] S. M. Berman, Phys. Rev. 112 (1958) 267. 

[85] T. Kinoshita and A. Sirlin, Phys. Rev. 113 (1959) 1652-1660. 

[86] G. Kallen, Springer Tracts in Modern Physics 46 (1968) 67-132. 

[87] T. van Ritbergen and R. G. Stuart, Phys. Rev. Lett. 82 (1999) 488, 
hep-ph/9808283 , 

[88] T. van Ritbergen and R. G. Stuart, Phys. Lett. B437 (1998) 201, 
hep-ph/9802341, 

[89] D. Bardin et al., "Z line shape", in Proc. of Workshop on Z Physics at LEP, 
Geneva, Switzerland, Feb 20-21 and May 8-9, 1989, report CERN 89-08 (1989) 
(G. Altarelli, R. Kleiss, and C. Verzegnassi, eds.), vol. 1, pp. 89-128. 

[90] D. Bardin, A. Leike, T. Riemann, and M. Sachwitz, Phys. Lett. B206 (1988) 
539-542. 

[91] F. A. Berends, W. L. van Neerven, and G. J. H. Burgers, Nucl. Phys. B297 
(1988) 429. 



59 



[92] B. Kniehl, M. Krawczyk, J. Kiihn, and R. G. Stuart, Phys. Lett. B209 (1988) 
337. 



[93] S. Jadach, M. Skrzypek, and M. Martinez, Phys. Lett. B280 (1992) 129-136. 

[94] D. Bardin, M. Bilenky, A. Sazonov, Y. Sedykh, T. Riemann, and M. Sachwitz, 
Phys. Lett. B255 (1991) 290-296. 

[95] M. Bilenky and A. Sazonov, "QED corrections at pole with realistic 
kinematical cuts", Dubna preprint JINR-E2-89-792 (1989). 

[96] P. C. Christova, M. Jack, and T. Riemann, Phys. Lett. B456 (1999) 264, 
|hep-ph/9902408 



[97] P. Christova, M. Jack, S. Riemann, and T. Riemann, "Predictions for fermion- 
pair production at LEP", in: J. Sola (ed.). Applications of Quantum Field 
Theory to Phenomenology, Proc. of 4th Int. Symp. on Radiative Corrections 
(RADCOR 98), Barcelona, Catalonia, Spain, 8-12 Sep 1998, p. 392-400 (World 



Scientific, Singapore, 1999), hep-ph/9812412 



[98] P. Christova, M. Jack, S. Riemann, and T. Riemann, "Predictions of ZFITTER 
v. 6 for fermion-pair production with acollinearity cut" , DESY preprint 99-037 



(1999),pIep^ph/9908289 



[99] D. Bardin, M. Bilenky, A. Chizhov, A. Sazonov, O. Fedorenko, T. Riemann, 
and M. Sachwitz, Nucl. Phys. B351 (1991) 1-48. 

[100] A. Leike, T. Riemann, and J. Rose, Phys. Lett. B273 (1991) 513-518, 
|hep-ph/9508390[ 

[101] T. Riemann, Phys. Lett. B293 (1992) 451-456, ple'p-ph/9506382| 

[102] M. Griinewald, S. Kirsch, and T. Riemann, Fortran program SMATASY version 
6.41 (18 Oct 2004), obtainable from 

/ afs / cern.ch / user /g/ gruenew / public / smatasy/smata6_41 .fort ran 
or from ^http : / / gruenew . home . cer n . ch /gruenew / [ 



[103] S. Kirsch and T. Riemann, Comp. Phys. Commun. 88 (1995) 89-108, 
|hep-ph/9408365[ 



60 



